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INTRODUCTION 


1.  LASER  TECHNOLOGY  AND  PROPAGATION 

The  computer  code  for  calculating  thermal  blooming  described  previously  has  been  extended 
to  include  the  effects  of  turbulence  in  the  atmosphere. 

A  theoretical  study  has  been  made  of  two  mechanisms  which  have  been  suggested  as  impor¬ 
tant  in  thermal  blooming  calculations:  the  dependence  of  a  upon  localized  temperature  in  a 
thermally  turbulent  medium,  and  the  effect  of  velocity  turbulence  which  would  tend  to  cool 
off  the  heated  air.  The  study  shows  that  for  power  levels  presently  contemplated,  both  these 
effects  are  negligible. 

A  serii  .<  of  propagation  calculations  have  been  performed  which  show  that  CVV  thermal  bloom¬ 
ing  can  lie  considerably  reduced  by  an  initial  phase  correction  of  the  laser  beam. 

Several  experiments  were  performed  to  investigate  the  effects  of  atmospheric  dust  on  the 
breakdown  threshold  of  10.6-p.ni  laser  beams  with  large  focal  spots.  The  lasers  used  in  these 
experiment::  were  a  1-kJ  pulsed  e-beam  at  Avco  Everett  Research  Laboratories  (AERL),  and 
a  3-|  pulsed  TEA  laser  at  Lincoln  Laboratory.  The  main  conclusions  from  these  experiments 
are:  the  breakdown  threshold  of  clean  air  is  close  to  that  predicted  by  microwave  breakdown 
theory;  dust  particles  can  lower  the  threshold  below  this  value  roughly  in  proportion  to  their 
size;  and  the  presence  of  atmospheric  dust  can  explain  the  "spot  size  effect"  reported  recently 
by  D.C.  Smith. 

Impulse  data  obtained  with  the  AERL  1-kJ  pulsed  CO.,  laser  have  been  roughly  correlated  with 
a  simple  cylindrical  blast  wave  model.  The  apparent  saturation  of  impulse  witli  increasing 
laser  energy  was  observed  in  the  presence  of  multiple  aerosol-induced  breakdowns  in  the 
laser  beam. 

Observations  of  plasmas  produced  by  irradiation  of  surfaces  with  the  AERL  pulsed  CO.,  laser 
were  carried  out.  The  plasma  and  associated  shock-wave  expansions  were  observed  with  a 

,  o 

framing  camera  (2  >  10  frames/sec).  The  plasma  absorption  at  wavelengths  4881  A  and 
10.6  pm  was  measured  and  used  to  infer  electron  density  histories. 

The  explosive  vaporization  of  individual  water  droplets  (5-  to  2 5 -pm  radius)  by  a  10.  6-pm 
laser  pulse  has  been  observed  with  a  high-speed  schlieren  photography  system.  The  hot  vapor 
sphere  and  the  shock  wave  produced  by  the  explosion  can  be  clearly  distinguished  in  the  schlie¬ 
ren  photographs. 

The  Lincoln  300-J  e-beam  excited  laser  is  in  the  final  stages  of  construction.  Laser  op¬ 
erations  should  commence  in  August  1972,  and  full  device  operation  is  expected  during  October. 

The  rate  equations  and  temperature  balance  equation  have  been  solved  to  predict  small  signal 
gain,  pulse  energy,  and  pulse  length  for  an  e-beam  excited  CO,,  laser.  The  initial  spike,  as 
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observed  experimentally  for  the  laser  pulse  from  the  Avco  e-beam  device,  has  been  predicted 
but  the  computed  temporal  width  is  narrower  th  .n  that  observed. 

A  2.5-m  unstable  confocal  positive  branch  resonator  has  been  designed  for  the  pulsed  500-J 
C02  laser  so  that  effective  use  can  be  made  of  the  large-volume  short-length  (1  m)  gain  me¬ 
dium.  Although  several  different  outcouplings  and  geometries  will  be  tried,  it  is  expected  that 
the  optimum  mode  quality  outcoupling  85  percent  will  produce  45  percent  of  the  net  energy 
within  the  first  Airy  dark  ring. 

The  Freedholm  integral  equation  is  used  to  determine  the  asymmetric  unstable  resonator  cav¬ 
ity  mode  field  distribution  without  a  gain  medium.  Conversion  is  made  by  function  expansion 
to  the  eigenvalue  solution  of  an  infinite  matrix.  A  finite  matrix  solution  approximation  is  then 
used  for  computer  calculation. 

A  3-nsec  600-keV  e-beam  source  has  been  used  to  initiate  a  CC>2  laser  beam.  Output  optical 
Pulses  up  to  120psec  and  22.  4J  have  been  observed  from  an  active  volume  of  about  1  liter. 

A  pulsed  C02  laser  beam  has  been  shaped  in  time  by  passing  the  beam  through  an  absorption 
cell  that  is  external  to  the  laser  and  filled  with  SF6  or  Freon-12.  The  initial  spike  on  the 
laser  pulse  has  been  successfully  attenuated  when  passed  through  the  absorption  cell. 

A  theoretical  investigation  has  been  carried  out  to  determine  the  difficulties  of  amplifying  very 
short  (~  Insec)  pulses  or  frequency-chirped  pulses  at  10-pm  optical  wavelength  to  high  ener¬ 
gies  (~200J/sequenec).  The  conclusions  indicated  that  only  the  volume  preionized  devices  are 
reasonable  in  the  last  stages  of  amplification.  Considerations  are  made  of  the  saturation  and 
pulse  sharpening  limitations  on  such  devices.  Frequency  distortion  and  index  changes  are  not 
seen  to  be  a  serious  problem.  The  most  severe  difficulties  are  found  to  be  in  hulk  medium 
heating  and  cathode  shock-wave  distortion  as  limiting  the  coherence  time  of  the  pulses, 

It  should  he  possible  to  produce  electron  densities  on  the  order  of  1012/em3  at  atmospheric 
pressure  by  photoionization  in  a  laser  gas  such  as  C02  or  CO.  The  ionization  produced  will  he 
uniform  throughout  large  laser  medium  volumes  if  the  photon  absorption  length  is  large,  and 
this  can  be  accomplished  by  seeding  the  laser  gas  with  a  low  partial  pressure  of  gas  that  has  a 
low  ionization  potential.  This  will  then  allow  deep  photon  penetration  and  low-energy  UV  pho¬ 
tons  to  be  utilized  for  ionization. 

A  small  Ionization  cell  containing  a  Xe  flash  lamp  with  suprasil  window  and  two  flat  parallel 
electrodes  is  used  to  study  volume  photoelectron  ionization.  When  0.  2torr  of  tripropylamine 
is  added  to  500-torr  N2>  an  electron  density  nQ  -  10*1  /cm3  is  measured  upon  flash-lamp  igni¬ 
tion.  Tlie  current  is  found  to  vary  linearly  with  applied  electrode  voltage,  which  implies  an  n 
production  only  by  photoionization. 

The  small  signal  gain  and  laser  level  pumping  rates  for  a  CC>2  laser  have  been  scaled  with 
pressureto  lOatm.  Using  longitudinal  excitation  by  a  104-A,  3-nsec  e-beam  source  (Ftbatron 
706),  the  peak  gain  for  the  P20  line  has  been  calculated  to  he  2.  2  percent/cm. 

3 

A  20-cm  C02  laser  device  has  teen  constructed  to  operate  at  lOatm,  and  experiments  with 
this  laser  will  commence  when  the  Felxitron  706  is  available  as  a  longitudinal  e-beam  source. 
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Measurements  and  evaluation  of  data  nave  been  completed  on  the  Lamb  dip  for  16  different 
rotational-vibrational  transitions  in  sealed-off  CO  lasers  operating  at  a  total  pressure  of 

1.5torr. 


II.  OPTICAL  MEASUREMENTS  AND  INSTRUMENTATION 

In  the  ZnS  upconversion  process,  AC  effects  have  been  measured  using  chopped  IR  radiation 
as  a  source  and  constant  UV  illumination  for  excitation.  The  AC  SNR  is  1  to  2  orders  of  mag¬ 
nitude  worse  than  the  SNR  obtained  in  previous  measurements  where  the  phosphor  is  exposed 
to  DC  IR  radiation  after  (not  during)  UV  excitation.  Apparatus  to  determine  the  S/N  charac¬ 
teristics  of  the  phosphor  screen  resolution  elements  is  being  constructed. 

Measurements  have  been  made  relevant  to  the  imaging  quality  of  a  cooled  4-inch-diameter, 
3-element  refractor  camera.  Tests  have  demonstrated  the  mechanical  integrity  of  the  design 
and  its  diffraction-limited  performance  in  the  wavelength  range  of  7.  5  to  12pm.  Extraneous 
radiation  at  angles  about  18  off  the  optical  axis  has  been  observed,  and  some  effort  has  been 
made  to  determine  the  cause  of  these  peaks. 

The  shearing  interferometer  was  employed  on  board  the  NKC-135  aircraft  with  auxiliary  op¬ 
tics  suspended  from  a  windew-mounted  air-foil  plate  to  measure  the  modulation  transfer  func¬ 
tion  of  the  local  boundary  layer  within  25cm  of  the  aircraft  skin.  The  existence  of  image  deg¬ 
radation-due  to  a  boundary  layer  was  established. 

Stellar  images  recorded  in  high-altitude  flight  with  a  diffraction-limited,  long-focal-length 
telescope  have  been  analyzed  to  determine  the  possible  effects  of  the  tropopause  in  image  deg¬ 
radation.  The  results  suggest  that  the  difference  between  the  stratosphere  and  troposphere 
is  not  measurable  in  the  presence  of  boundary-layer  turbulence,  and  a  determination  of  turbu¬ 
lence  effects  agrees  quite  well  with  interferometer  results  measured  separately. 

Detailed  tunable  diode  laser  studies  have  been  made  of  the  j-  band  of  SO.,  centered  around 
8.7pm.  Comparison  between  the  experimental  results  and  theoretical  predictions  based  on 
microwave  data  has  resulted  in  the  establishment  of  a  more  accurate  value  for  the  band  center 
frequency.  Characteristics  of  the  absorption  spectra  under  varying  conditions  of  temperature 
and  pressure  arc  used  to  determine  optimum  wavelengths  for  monitoring  applications. 

Several  extensions  of  tunable  diode  laser  techniques  to  air-pollution  monitoring  are  described. 
The  sensitivity  of  acousto-optic  absorption  cells  and  direct  transmission  measurements  for 
air-pollution  monitoring  have  been  compared;  real-time,  quantitative,  in-situ  measurements 
of  auto  exhaust  have  been  made  and  are  reported;  and  long  path  laser  transmission  studies  of 
ambient  air  pollution  are  described. 


III.  LASER  IMAGING  RADAR 

The  10. 6-pm  imaging  radar  has  been  used  to  measure  cross-section  characteristics  of  two 
satellite  models  and  a  Terrier  radome. 
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A  simple,  portable  modulation-transfer-function  (MTF)  Instrument  based  on  a  Sieman’s  star 
reticle  has  been  built  and  tested.  Under  laboratory  conditions,  the  instrument  has  demon¬ 
strated  diffraction-limited  performance,  in  the  field,  the  atmospheric  MTF  has  been  meas¬ 
ured  under  a  wide  variety  of  conditions  in  support  of  atmospheric  spatiai  coherence  measure¬ 
ments  which  were  being  made  with  a  48-inch  telescope  at  10.6pm. 

Electronics  have  been  constructed  both  for  control  of  the  MTF  instrument  and  for  processing 
of  the  output  signal.  The  MTF  processor  provides  a  calibrated  output  voltage  proportional  to 
the  percentage  modulation.  Various  control  circuits  for  the  Firepond  laser  radar  have  also 
been  designed  and  implemented. 

Spatial  coherence  measurements  have  been  made  at  10.  6  and  0.  63  pm.  The  data  at  10.6pm 
are  currently  being  reduced,  but  qualitatively  agree  with  the  results  obtained  at  0.63pm. 
These  latter  data,  taken  over  a  5.4-km  horizontal  path  at  an  average  68m  above  the  ground, 
indicate  that  the  coherence  diameter  rQ  (0.63pm)  for  a  wide  variety  of  atmospheric  conditions 
varies  from  1cm  <rQ  (0.63pm)  <  20cm.  From  this,  we  predict  that  the  10.  6-pm  data  will 
reveal  30cm  <  rQ(10.6pm)  <  6  m  for  this  path. 

A  10.  6-pm  isolator  has  been  fabricated  using  Faraday  free-carrier  rotation.  The  isolator  has 
a  forward  insertion  loss  of  2dFi,  a  total  reverse  loss  of  22 dB,  and  a  maximum  throughput 
power  of  30  W  CW. 
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I.  LASER  TECHNOLOGY  AND  PROPAGATION 


A.  CW  THEORY 

1.  Laser  Propagation  Through  a  Turbulent  Medium 

In  previous  reports,  we  discussed  the  propagation  of  a  high-power  laser  beam  through  a 
nonturbulent  atmosphere.  In  the  past  six  months,  we  have  extended  our  capabilities  to  include 
propagation  through  a  turbulent  atmosphere. 

The  standard  approach  to  the  calculation  of  propagation  of  a  low-power  beam  through  a  tur¬ 
bulent  medium  is  by  analytic  computation  of  ensemble  average  values  of  various  quantities.  (A 
good  review  of  this  field  is  found  in  Ref.  1.)  However,  since  analytic  methods  had  not  been  suc¬ 
cessful  in  the  high-power,  nonturbulent  case,  another  approach  was  needed.  We  accomplished 
this  by  simulating  individual  realizations  of  the  atmosphere,  with  the  constraint  that  the  ensem¬ 
ble  of  these  realizations  should  have  the  statistical  properties  of  an  accepted  (or  at  least,  widely 
used)  model  of  the  turbulent  atmosphere.  In  the  following  section,  we  argue  that  the  velocity 
turbulence  of  the  atmosphere  is  of  little  importance  except  in  the  case  of  an  extremely  high- 
power  beam,  and  that  for  the  majority  of  cases  it  is  an  adequate  approximation  simply  to  add  the 
naturally  occurring  temperature  (or  index)  fluctuations  to  the  temperature  (or  index)  change  in¬ 
duced  by  the  beam  itself.  Since  our  numerical  method  is  capable  of  computing  the  propagation 
through  an  arbitrary  atmosphere  (so  long  as  the  index  gradients  are  not  excessive),  the  only 
remaining  problem  is  the  generation  of  suitable  realizations  of  the  temperature  fluctuations. 

For  our  purposes,  sufficient  information  about  the  statistical  properties  of  the  index  of  re¬ 
fraction  is  contained  in  the  autocorrelation  function  <n(  r)  n(  r  +  p  )>.  If  the  atmosphere  is  con- 
sidered  to  be  statistically  homogeneous,  this  function  is  independent  of  r,  and  can  be  described 
equally  well  by  the  (three-dimensional)  spectrum 


4>,(iT)  =  — L r  f  dVeik'P<n(r)n(r+p)> 

J  (2*)  J 

This  spectrum  is  generally  used  for  the  statistical  description  of  the  atmosphere.  But  in  our 
propagation  calculations,  the  quantity  n  does  not  enter  directly,  but  only  as  an  integral,  f  ndz. 
To  this  integral  corresponds  a  two-dimensional  spectrum: 


$2fkx’  ky*  =  j-^2  ^  dxdy  exp[i(kxx  +  ky>')]  1  n(r)dzj  1  n(r+p’)dz^  . 


.  7, .  +Az 


It  is  easy  to  show  that  $2  *s  related  to  $  ^  by  the  equation 

$,(k  ,  k  )  =  (A z  )^  \  dk  <J>  7(k  ,  k  ,  k  )  fj  (k  Az/2)]^ 
2  x  y  J  z  3  x  y1  z  lJo  z  '  ‘ 


where  jQ(u)  =  sinu/u. 
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Fig.  1-1 .  Iso-irradiance  contour  plots  for  conditions  corresponding  to  Firepond  focility:  122-cm 
Cassegroinion  (0. 33  obscuration  ratio)  focused  at  5.4  km,  \  =  10.6  pm.  Turbulence  is  assumed 

to  be  homogeneou'  and  isotropic,  with  =  2  X  10  Four  different  realizations.  Grid 

spocings  ore  10  cm. 
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We  now  wish  to  generate  values  of  j ndz  on  a  finite  two-dimensional  grid.  We  do  this  by 

generating  an  array  of  random  complex  numbers  qm  n  with  zero  mean  and  Gaussian  distribution 

in  the  space  k  k  .  multiplying  by  the  required  function  of  k  ,  k  ,  and  taking  a  Fourier  transform: 
x  j  x  y 


ndx 


(2nf  v 


L  L 
x  y 


exp 


m,n 


L  .  / mx 

„  /27rm 
32\  L  ' 

27rn\ 

L  I 

1  '  X 

y/J  Sj 

'  X 

y' 

where  and  L  are  the  total  widths  of  the  grid  in  the  x-  and  y-directions,  respectively. 

Requiring  a  least-mean-square  fit  of  J  ndz  based  on  a  finite  grid  to  the  continuum  expres¬ 
sion,  one  can  show  that  a  convolution  is  necessary  in  order  to  take  proper  account  of  frequencies 
present  in  the  continuous  spectrum  but  not  present  in  the  spectrum  S.,  for  the  discrete  wave- 

numbers  k  =  2jrm/L  and  k  =  2irn/L  .  S?  is  then  related  to  by 
a.  x  y  y  c. 


S,(k  ,  k  ) 
2  x  y 


L2  L2 

X 


(27 r 


MI 


dk'  dk1  a  <k'  ,  k'  ) 
x  y  2  x*  y' 


(k.. 


k1  )  L 
x  x 


(k  -  k'  )  L 

y  y  y 


We  use  the  scheme  described  above  to  generate  an  array  of  phase  changes  (proportional  to 
f  n a/)  twice  for  every  step  in  z,  in  accordance  with  our  method  of  calculation.2  The  arrays  are 
taken  to  be  statistically  independent  of  one  another,  since  the  z-steps  actually  used  are  consider¬ 
ably  larger  than  the  outer  scale  of  turbulence,  except  in  the  immediate  vicinity  of  a  focus  where 
the  small  optical  lever  arm  makes  phase  changes  of  little  importance.  For  the  same  reason, 
we  believe  that  the  assumption  of  a  Gaussian  distribution  for  the  phase  changes  is  justified.  For 
the  spectrum  <k3,  we  have  used  the  modified  von  Karman  spectrum  in  the  form 

*3<k)  =  0.033  C2  [k2  +(^l)2]‘11/6  exp (-  (kfQ/2*)2] 

where  Lq  is  the  outer  scale,  and  is  the  inner  scale  of  turbulence. 

Figure  1-1  shows  some  iso-irradiance  contour  plots  that  result  when  we  apply  our  method 
to  a  particular  experimental  situation,  without  nonlinear  effects.  Each  plot  corresponds  to  an 
instantaneous  picture  of  the  beam,  and  the  variation  from  one  plot  to  another  is  due  to  the 
statistical  fluctuations  of  the  assumed  atmosphere. 

As  overall  measures  of  the  effect  of  turbulence,  we  have  chosen  two  quantities.  One  is  the 
equivalent  area,  defined  as 

A  _  (  f  ldxdy)2 
eq  2  /  I2  dxdy 

For  well-behaved  intensity  patterns,  the  equivalent  area  so  defined  agrees  with  the  area  defined  by 

Arms  =  *  ^  [<x  “  x  )2  +  (y  -  y  )2]  ldxdy/  (  ldxdy 

based  on  the  standard  deviations.  But  for  the  Airy  pattern,  the  rms  area  diverges  whereas  the 
equivalent  area  gives  a  reasonable  value.  In  addition,  the  standard  deviations  cannot  be  reliably 
calculated  for  noisy  experimental  data,  since  the  weignting  of  the  intensity  with  r3  (as  can  be 
seen  by  putting  the  integral  into  cylindrical  coordinates)  leads  easily  to  numerical  difficulties 
for  large  r.  We  therefore  prefer  to  use  the  equivalent  area.  This  is  usually  divided  by  the 
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corresponding  quantity  in  the  absenee  of  turbulence  Aq,  and  a  nondimensional  "beam  spread  is 
defined: 


a?  =  (A 
A  eq 


A  )/A 
o  o 


The  other  quantity  is  the  displacement  of  the  eenter  of  irradiance  of  the  pattern: 


=  xldxdy/^  Idxdy 


w 

y 


1  dx  dy 


and  the  nondimensional  "wander"  (squared)  is  defined  as 


ffw  E*(wx2  +wy)/Ao  • 

2 

In  Figs.  1-2  and  1-3,  we  plot  the  values  of  these  quantities  as  a  function  of  Cn>  the  measure  of 
the  strength  of  turbulenee.  At  a  given  value  of  C'n ,  the  different  points  correspond  to  different 
realizations  of  the  atmosphere.  Points  that  are  eonneeted  by  dashed  lines  correspond  to  reali¬ 
zations  that  differ  only  by  the  value  of  C^;  i.e.,  the  random- numbers  used  to  generate  the  reali¬ 
zations  are  the  same. 

The  fact  that  these  dashed  lines  are  nearly  parallel  to  one  another  suggests  that  it  maybe 
possible  to  study  the  dependence  of  turbulent  propagation  on  various  parameters  by  maintaining 
a  single  set  of  random  numbers  and  varying  the  parameters,  without  having  to  accumulate  a 
large  number  of  statistics  for  each  set  of  parameters.  This  possibility  is  made  more  plausible 
by  an  examination  of  the  eontour  plots  of  Figs.  1-4  and  1-5,  whieh  result  from  just  such  a  study, 
where  a  set  of  dimensionless  numbers  appropriate  for  this  problem  have  been  chosen  and  varied 
systematically.  The  dimensionless  numbers  are  the  Fresr.el  number 


ttie  f/number 

f/  =  lV2ani 

the  "turbulent  distortion  number" 


NT  =Vk'CnljV5 


and  the  ratios  a  /L  and  a  fl  .  Here,  a  is  a  eharaeteristie  radius  in  the  initial  plane,  k  is 
m  o  rn  o  ni  ,  «  » 

the  optical  wavenumber,  and  both  the  initial  radius  of  eurvature  and  the  distance  to  the  plane  of 

observation  are  taken  to  be  U  For  Figs.  T-4  and  1-5,  the  parameter  N,f  has  been  kept  constant, 

while  the  other  parameters  are  varied.  From  the  "family  resemblance"  of  the  plots,  it  seems 

clear  that  scaling  laws  ean  be  deduced  in  this  way  and  that  N,j,  is  the  most  important  parameter 

in  the  range  studied.  Its  effect  on  the  beam  spread  and  wander  are  already  visible  in  Figs.  1-2 

and  1-3;  Figs.  1-6  and  1-7  show  the  effeets  of  varying  Nj,. 

We  can  now  proeeed  to  eompute  the  propagation  of  a  high- power  beam  through  a  turhulent 
atmosphere.  In  Fig.  1-8,  we  see  a  comparison  of  results  for  a  high-power  beam  propagating 
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Fig.  1-4.  Iso-iriodionce  contours  for  geometry  of  Fig.  1-1 ,  with  NT  =  1 . 86  ond  vorious  volues 
of  Nc  ond  f/.  1 
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Nf  =  2  0 . 4  f/6270 


Fig.  1-5.  Iso-irrodionce  contours  for  geometry  of  Fig.  1-1 ,  with  NT  =  3.  83  ond  vorious  volues 
of  Np  end  f/. 
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Fig.  1-8.  Thermol  blooming  in  presence  of  turbulence.  Top  contour  plot  shows  beom  with  turbulence 
absent  ond  thermal  blooming  present.  Left  three  contour  plots  show  beam  in  obsence  of  thermol  bloom¬ 
ing,  ofter  propagation  through  three  realizations  of  turbulent  atmosphere  (C2  =  10"'4  m_2/3).  Right 
three  contour  plots  show  thermally  blooming  beam  ofter  possing  through  same  three  realizations  of  tur¬ 
bulent  otmosphere.  Initial  beom  was  on  infinite  Gaussion,  l/e2  power  dicmeter  1  m,  focused  ond  ob¬ 
served  at  ronge  of  10  km,  X  =  10.6  pm,  ond  oP/v^  ~  10"4. 
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Fig.  1-9.  Total  loss  (thermol  blooming  plus  turbulence)  as  a  function  of  slewing  rate.  Solid  line  shows 
loss  in  absence  of  turbulence,  while  hotched  oreas  show  range  of  loss  ossocioted  with  various  values  of 
.  Infinite  Goussian,  1/e2  power  diometer  1  m,  ronge  3  km,  X  =  20  pm,  oP/v  =  10-4. 


through  the  same  turbulent  atmosphere.  Figure  1-9  shows  the 
and  thermal  blooming  are  combined,  for  a  particular  example, 
under  way. 


intensity  loss  when  turbulence 
More  detailed  studies  are 

L.  C.  Bradley 
J.  Herrmann 


2.  Interaction  of  Heating  and  Turbulence  in  Thermal  Blooming 

In  our  calculation  of  the  propagation  of  a  high-power  beam  through  a  turbulent  atmosphere, 
we  assumed  that  the  temperature  variations  due  to  heating  may  be  simply  added  to  those  already 
present  in  the  atmosphere.  It  has  been  suggested  that  the  interaction  between  heating  and  tur¬ 
bulence  may  give  rise  to  additional  temperature  variations  sufficiently  large  that  the  propagation 
will  be  significantly  affected.  Two  mechanisms  have  been  proposed: 

(a)  Warmer  blobs  have  a  larger  absorption  coefficient  than  cooler  ones, 
so  heating  will  tend  to  amplify  the  variations. 

(b)  Velocity  fluctuations  in  the  medium  will  tend  to  cool  off  the  heated 

air  more  effectively,  and  also  will  cause  further  temperature  fluctuations. 

Our  purpose  here  is  to  estimate  these  effects  and  show  that  both  are  negligible,  at  least  at  pres¬ 
ently  contemplated  power  levels. 

We  assume  that  wind  velocities  arc  much  smaller  than  the  sound  velocity,  so  that  the  pres¬ 
sure  may  be  considered  constant.  The  temperature  is  thus  considered  a  passive  additive  to  the 
turbulent  velocity  field.  If  the  winds  arc  not  so  slow,  the  heating  effects  will  be  negligible  any¬ 
how,  Then,  the  equation  of  energy  conservation  may  be  written  as 
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%£  +  v  .  VT  =  — 
9t 


(1-1) 


We  assume  the  properties  of  the  radiation  source  to  be  independent  of  time  (or  slowly  changing), 
and  divide  the  properties  of  the  medium  into  mean  and  fluctuating  parts: 


T  -  x  +  T. 
o  1 


V  -  V  +  V, 

o  1 


-  = (-)  + (-) 
p  Vo  +  Vi 


The  mean  parts  are  assumed  to  satisfy 


v  .  VT  =  (2.)  — 
o  o  p  o  c 


(1-2) 


To  this  might  be  added  a  term  -v.  •  VT^,  but  we  argue  below  that  this  term  is  negligibly  small 
for  practical  application.  Then  we  have 


3T. 


~  +  v  •  VT  +  v,  •  VT,  =  (-),  —  -  v, 
St  o  t  1  1  p  1  c  1 

p 


VT  - 


ST 

c 

St 


d-3) 


We  now  separate  that  part  of  the  random  temperature  field  which  can  be  considered  as  the  frozen 
turbulence  T^.,  satisfying  the  Taylor  hypothesis 

ST  _ 

It"  ="V0  •  VTf  • 


Then,  with  -  Tj.  +  we  get 
ST 


vp  +  (  v  +  v,  )  •  VT,  =  ( - )  —  -  v.  •  VT  -  v, 

St  o  1  2  p  1  c  1  o  1 

P 


ST 

VT - — 

v  f  St 


d-4) 


In  our  propagation  code  (see  Sec.  A-l )  in  a  turbulent  atmosphere  the  contributions  of  T  and  T„ 

O  i 

have  been  included.  We  must  now  estimate  the  contributions  on  the  right-hand  side  of  Kq.  (1-4) 
to  show  that  T^  is  negligible. 

The  time  for  a  fluid  particle  to  cross  the  beam  is  ~a/v  ,  where  a  is  the  beam  radius,  so 
the  increment  in  T^  is 


a^ 

2  ~  v. 


-)  - 
pic 


ST 
_ r 

St 


We  may  compare  the  increment  from  the  first  term  (AjT^  to  the  rms  value  of  Tf  in  the  unper¬ 


turbed  fluid,  which  we  call  T  : 

’  r 


Now, 


A  T 

1  2  J_  a_  .  a,  J_ 

T  ~  T  v  '  p’l  c 
r  r  o  p 


it!l  „  J1  Zr  „  6  ,  Ll 

,  (V.  kT  T  D'  ~  T 

(p>o 
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inhere  E  is  the  lower  level  of  the  laser  transition,  and  I  -  P/Va  ,  so  we  find 


AfT2 


6.  5cv  P 
_ o 

7rp  c  T  v  a 
o  p  o  o 


For  fairly  extreme  cases,  this  may  amount  to  percent,  which  we  consider  to  be  negligible. 

It  is  not  appropriate  to  compare  the  increment  from  the  second  term  (A^T^)  with  T  ,  since 
under  conditions  of  low  lapse  rate  (at  dawn  or  dusk)  it  is  possible  to  have  large  velocity  fluctua¬ 
tions  and  small  temperature  fluctuations.  Instead,  we  compare  it  with  AT^,  the  change  in  Tq 

in  traversing  the  beam.  It  is  easily  seen  that  (A,T,/AT  )  ~  (v  /v  ),  where  v  is  the  rms  value 
h  J  2  2'  o  r  o  r 

of  v,.  This  is  an  overestimate  of  the  effect,  however,  for  the  important  quantity  for  propagation 

is  not  the  change  in  temperature  itself,  but  its  integral  over  the  length  of  the  beam,  which  is 

proportional  to  the  phase  change.  Since  v^  is  uncorrelated  over  paths  greater  than  some  outer 

scale  Lq,  while  adds  coherently  over  the  whole  path,  which  we  assume  is  much  larger  than 

L  ,  we  have 
o’ 

(  AT^dz  oc  (ATo)  (Az) 


A2T2dz=c(A2T2>(I-oA-> 


1/2 


and  the  ratio  of  the  corresponding  phase  changes  is 

V 

A<po  ~  vo  lA'xJ 

At  most,  (v  /v  )  ~  0.1,  while  (I.^/Az)  ~  0.01  for  the  cases  in  which  we  are  interested.  Thus, 
wo  consider  that  this  term  too  can  he  neglected. 

The  third  term  (A^T^)  describes  the  breakup  of  frozen  turbulence.  Since  VTj,  ~  (T^/l,  ), 
we  find 


A3T2 


v 

r 


a 

L 

o 


which  is  negligibly  small. 

The  fourth  term  (A.T_)  is  the  most  difficult  with  which  to  deal.  The  time  variation  of  T  is 
4  2  o 

due  to  the  fact  that  turbulence  in  the  initial  parts  of  the  beam  leads  to  scintillations  (time 
changes  in  I)  at  later  points  in  the  beam.  In  general,  the  rapidly  varying  part  of  I  has  high 
spatial  frequencies  and  therefore  tends  to  integrate  out  in  the  equation  for  T^,  whereas  the  low 
spatial  frequeheies,  which  are  the  most  important  in  contributing  to  T^,  are  slowly  varying. 
Furthermore,  in  a  focused  beam,  the  relevant  time-scale  a/v^  becomes  smaller  as  the  beam 
travels  farther  and  as  the  scintillations  thus  become  larger.  Therefore,  A^T^  may  not  be  as 
important  as  it  appears  at  first  glance.  Tn  an  attempt  to  make  these  statements  more  quantitative, 
we  are  now  carrying  out  "computer  experiments"  on  low- power  beams. 

The  above  considerations  apply  to  a  beam  of  moderate  dimensions  (<f  m)  propagating  through 

3 

the  atmosphere.  Recently,  Gebhardt,  £t  ah  have  studied  thermal  blooming  in  a  turbulent  me¬ 
dium  in  the  laboratory.  They  have  attempted  to  make  v  and  T,  as  small  as  possible,  and  the 

o  !  ^ 

integral  scale  for  their  turbulent  medium  is  on  the  order  of  a  beam  diameter.  For  such  an 
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environment,  it  is  no  longer  permissible  to  neglect  the  term  -v^  •  VT^  in  Eq.  (1-2);  the  usual 
practice  is  to  replace  it  by  an  eddy  diffusion  term  V  •  (KeVTQ)  (see  Ref.  5).  While  this  is  no 
doubt  appropriate  for  their  experiment,  we  believe  that  our  order-of-magnitude  analysis  above 
shows  that  this  term  will  not  be  important  in  the  open  atmosphere. 

J.  Herrmann 
L.  C.  Bradley 

3.  Calculations  of  Phase  Correction  for  a  Laboratory  Thermal 
Blooming  Experiment 

In  order  to  design  n  laboratory  experiment  for  the  study  of  the  elimination  of  thermal  bloom¬ 
ing  by  an  initial  phase  correction,  a  series  of  propagation  calculations  has  been  performed. 

Figure  I- 1 0  shows  the  radius  of  the  unbloomcd  beam  near  t'  focus  for  six  aperture  sizes. 
All  cases  assume  a  Gaussian  beam  truncated  at  this  aperture  si  with  a  radius  of  curvature  of 
1  m.  From  Fig.  1-10,  the  optimum  experimental  condition  maybe  determined. 


Fig.  1-10.  Equivolent  rodius  of  unbloomed  beom 
for  s'x  operture  sizes.  Goussion  beam  is  truncated 
at  l/e^  power  level  and  has  wavefront  curvoture 
of  1-m  rodius.  All  curves  apply  to  wovenumber 
k  =  1.2  /  105  cm-1. 


Figure  1-11  shows  the  5- mm-aperturc  beam  with  thermal  blooming  (curve  B).  By  applying 
an  initial  phase  correction,  thermal  blooming  is  considerably  reduced  (curve  ('),  approaching 
the  unbloomcd  cond  on  (curve  A). 

Figure  1-12  shows  the  equivalent  curves  fora  2-mm-aperture  beam. 

J.  Herrmann 
1„  C.  Bradley 
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Fig.  1-11.  For  a  beam  with  a  5-mm-diameter 
aperture,  curve  A  is  the  unbloomed  case,  curve  B 
is  with  thermal  blooming,  and  curve  C  is  with 
initial  phase  correction.  Blooming  parameters 
are:  P  =  1  W,  a  =  10"  3  cm-' ,  v  =  0.2  cm/sec, 
and  u  =  0.2  rod/sec. 
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Fig.  1-12.  Same  as  Fig.  1-1 1  for  2-mm-aperture 


beam.  Blooming  parameters  a!sa  same,  exce 
v  =  0.5  cm/s ec  and  u  =  0.5  rad/sec. 


Fig.  1-13.  Breakdown  thresholds  measured  for  spot 
diameters  0. 3  and  0.5  mm.  Data  for  smaller  spot 
size  are  from  TEA  laser  results  at  Lincoln  Labora¬ 
tory  and  from  Ref.  10. 
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B.  EFFECTS 

1.  Laser-Induced  Air  Breakdown 

Several  experiments  were  performed  to  investigate  the  effects  of  atmospheric  dust  on  the 
breakdown  threshold  of  10.6-pm  laser  beams  with  large  focal  spots. 

a.  A  ERL  Breakdown  Experiment8,  7 

In  March  1972,  we  conducted  a  series  of  experiments  at  Avco  Everett  Research  Laboratories 
(AERL)  using  their  1-k.J  pulsed  e-beam  laser.  This  laser  emitted  a  10-  x  20-em  beam  which 
was  from  1  to  2X  diffraction-limited  and  had  a  pea«  pulse  power  »4  x  108W.  With  ?5-m  focusing 
optics,  we  were  able  to  obtain  breakdown  at  a  maximum  spot  area  of  »0.  3  cm  .  The  focused 
power  was  varied  by  either  changing  the  sustainer  voltage  or  by  inserting  finely  perforated 
screens  with  geometrical  transmission  of  0.70  and  0.5  in  the  beam.  The  threshold  for  breakdown 
was  obtained  to  within  a  factor  of  approximately  two.  The  results  for  Inrge  spot  sizes  are  shown 
in  Eig.  1-13  together  with  Lincoln  laboratory  TEA  laser  results  and  the  results  from  D.  C.  Smith8 
As  can  be  seen,  the  threshold  continues  to  decrease  at  these  spot  sizes. 

In  general,  for  flux  densities  larger  than  the  threshold,  breakdown  was  not  confined  to  the 
focal  volume  but  instead  consisted  of  a  string  of  discrete  breakdowns  extending  up  the  beam 
(Eig.  1-14).  The  length  of  this  string  increased  with  increasing  laser  power,  as  shown  in  Fig.  1-15. 

The  breakdown  string  always  extended  well  beyond  the  point  where  the  beam  flux  was  equal  to  the 

9  2  ^ 

clean-air  breakdown  value  («3  x  10  W/cm  ).  The  minimum  flux  densities  where  breakdown  oc¬ 
curred  were  «io7W/cmt. 

These  results  can  be  understood  by  assuming  the  discrete  breakdowns  are  ocmirring  off 
large  dust  particles  in  the  beam.  The  size  distribution  of  dust  in  the  experimental  area  was 
monitored  using  a  light  scattering  particle  counter  (Royco  Model  225)  and  a  256-channel  pulse- 
height  analyzer  (Northern  Model  NS-6  33).  Particles  in  the  size  range  1  to  12  pm  could  be  meas¬ 
ured  in  this  way.  A  typical  measured  distribution  is  shown  in  Fig.  1-16. 


Fig.  1-14.  String  of  breokdowns  extending 
from  focol  spot  back  toward  laser. 
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Fig.  1-17.  Beom  geometry  ond  nototion.  Breokdown  string  extends  to  position  Z 
where  flux  density  is  I(Zq). 

The  functional  dependence  of  the  breakdown  threshold  on  particle  size  can  be  obtained  in  the 
following  way:  f^t  IfZ^)  be  the  flux  density  at  the  furthest  extent  of  the  breakdown  string  Z^ 

(Fig.  1-17).  Assume  that  the  breakdown  is  determined  by  the  largest  particle  size  D  that  can  be 
present  near  this  position.  The  density  of  particles  of  size  or  larger  is  given  by 

MDo)=!”  ^dDarlO-^D^cm-3  . 
o 

The  largest  partiele  that  can  be  present  near  Zq  is  given  by  the  conditions 
N(I)0)  V(Zo)  >10 

where  V(Z  )  is  the  volume  of  the  laser  beam  from  the  focal  spot  to  Z  .  That  is.  there  will  be  a 
o  o 

particle  of  size  Dq  near  7.Q  if  there  are  approximately  ten  such  particles  in  V(Zq),  This  gives 
l)Q  -  f(ZQ)  which,  together  with  IfZ^),  gives  the  threshold  intensity  on  partiele  size  shown  in 
Fig.  1-18.  Note  that  the  results  for  I)  are  not  very  sensitive  to  the  number  of  particles  assumed 
to  be  present  in  VfZ^). 

We  conclude  that  the  experimental  results  for  breakdown  threshold  can  be  understood  if  the 
breakdown  threshold  depends  on  aerosol  particle  size  in  the  manner  shown  in  Fig.  1-18. 


PARTICLE  DIAMETER  D  lum) 
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b. 


Breakdown  Results  with  TEA  leasers 


7 


It  has  been  reported  recently  that  the  breakdown  threshold  for  10.6-pm  radiation  decreases 

with  increasing  spot  size  well  below  the  value  predieted  by  classical  microwave  breakdown 

-2  -1 

theory.  We  have  verified  the  effeet  for  focal  spot  sizes  in  the  range  10  to  1  0  cm,  but  found 
that  the  threshold  was  sensitive  to  the  ;•  nount  of  dust  near  the  focal  volume,  especially  at  large 
spo*  sizes. 

A  chamber  was  designed  and  constructed  which  provided  a  variable  atmosphere  from  dusty 
to  highly  filtered  clean  air.  The  dust  count  and  size  distribution  within  the  chamber  were  mon¬ 
itored  with  a  light  scattering  particle  eounter  (Royco  225)  and  a  256-channel  pulse-height  analyzer 
(Northern  NS-633).  We  were  thus  able  to  measure  the  size  distribution  in  the  range  1  to  12  pm. 

A  series  of  threshold  measurements  was  conducted  in  known  atmospheres.  For  the  very- 
clean-air  case,  the  threshold  approached  an  asymptotic  value  at  large  spot  size  which  was  close 
to  that  expected  from  microwave  breakdown  theory.  Figure  1-19  shows  these  data  along  with  a 
theoretical  curve  obtained  from  Ref.  9.  The  focused  spot  was  reetangular  with  dimensions 
S  x  5S.  The  fit  to  the  data  was  for  a  diffusion  length  A  =  S/6.  The  other  data  in  Fig.  1-19  corre¬ 
spond  to  an  atmosphere  of  carbon  dust.  As  can  be  seen,  the  threshold  continues  to  decrease 
with  increasing  spot  size.  The  relative  dust  conditions  for  the  three  cases  are  indicated  in  the 
figure  by  the  total  density  of  partieles  with  diameters  larger  than  2  pm. 

As  in  the  previous  section,  these  results  can  be  explained  by  assuming  that  the  breakdown 
threshold  is  determined  by  the  largest  particle  present  in  the  foeal  volume.  Thus,  for  example, 
at  S  =  0.03  cm  in  Fig.  1-20  the  three  dust  conditions  give  maximum  partiele  sizes  of  3,  1,  and 
«1  pm  for  the  three  data  points.  The  inferred  threshold  dependence  on  particle  size  obtained  in 

this  way  is  shown  in  region  A  of  the  figure.  Also  shown  are  the  results  discussed  in  the  previous 

1 0 

section  and  the  thresholds  measured  directly  for  thin  glass  fibers  and  single  particles. 

The  above  results  indicate  that  atmospheric  aerosols  can  aeeount  for  the  threshold  lowering 
at  large  spot,  sizes.  Preliminary  calculations  of  the  cascade  proeess  in  exploding  particle  vapors 
give  the  right  order  of  magnitude  as  well  as  the  particle- size  dependence  of  the  breakdown 
threshold. 

Further  work  in  this  area  will  inelude  studies  of  different  material  aerosols,  some  refine¬ 
ment  of  the  theoretical  models,  and  studies  of  the  interaction  of  radiation  with  aerosols  below 


the  breakdown  threshold. 


D.  E.  Lencioni 
L.  C.  Marquet 


2.  Surface  Interactions 

a.  Impulse  Delivery  by  a  Pulsed  Laser  at  10.6  pm 

Measurements  and  a  correlation  of  impulse  delivered  to  targets  of  various  sizes  and  ma¬ 
terials  with  the  use  of  the  AERL  1-kJ  CO_  laser  are  reported  in  this  section.  Some  of  these 

|  j 

measurements  and  the  experimental  procedure  are  reported  in  detail  elsewhere. 

The  total  impulse  delivered  to  carbon  a  id  aluminum  targets  of  different  sizes  was  measured 

2  12 

for  focal  spot  areas  of  approximately  0.1  and  0.3cm  (27-m  focal-length  mirror).  AERL  data 

2 

for  aluminum  and  tungsten  targets  are  also  reported  for  focal  spot  areas  of  0.02  and  0.6  cm 
(10.6-m  focal-length  mirror).  The  delivered  impulse  was  found  to  be  independent  of  the  target 
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material  (carbon,  aluminum,  and  tungsten).  There  was  no  measurable  mass  loss  from  the 

-4  -  3 

aluminum  targets  (Am  <10  g)  and  the  mass  loss  from  the  earbon  targets  was  less  than  10  g. 

2 

The  projected  target  areas  ranged  from  0.7  to  40cm  . 

A  simple  model  based  on  an  expanding  cylindrical  blast  wave  and  a  plasma  front  velocity 
proportional  to  the  eube  root  of  the  power  density  is  used  to  obtain  impulse  sealing  relations.  It 
is  shown  in  Refs.  13  and  14  that  the  eylindr ieal-blast-wave  relations  ean  be  modified  to  approx¬ 
imately  account,  for  time-dependent  energy  deposition.  The  shock-wave  radius  and  the  resulting 
pressure  are  shown  to  be  proportional  to  the  following  relations: 


R(t)  « 


(1-5) 


P(t)  «  K|(t) 


-1 


(1-6) 


where  the  energy-per-unit-length  K.  is  defined  as 

F  K(t) 

f  "  nt)  ■ 


(1-7) 


Here,  K(t)  is  the  total  energy  deposited  in  the  plasma  and  f(t)  in  the  distance  the  plasma  front  has 
moved  away  from  the  target  during  the  time  t.  Then,  for  the  ease  of  a  eonstant  power  density, 
we  find 


“  Ao(q)2/3 


(1-8) 


where  is  the  foeal  spot  area  and  q  is  the  average  power  density  incident  on  the  target.  The 
impulse  is  given  by 


1  «  P(t)  R(t)2  dt 
'  o 

and  for  the  case  of  a  eonstant  power  density 
IocAo(q)2/3r  . 


(1-9) 


(1-10) 


Two  limiting  eases  may  be  identified  which  determine  the  termination  of  the  impulse  delivery  r. 
For  small  free-standing  targets,  the  impulse  delivery  ends  when  the  cylindrical  shock  wave  ex¬ 
pands  off  the  target.  In  this  ease,  Fq.  (1-5)  yields  r  re  A,rA^  /2(q  )_1  '3,  where  AT  is  the  target 
area.  The  impulse  sealing  relation  is 

1  “  ATA01//2(q)1//3  •  (Ml) 


For  an  infinite  target,  the  impulse  delivery  ends  when  the  shoek-wave  pressure  has  decreased 
to  about  1  atm.  Then,  Fq.  (1-6)  yields  r  «  A^^fq)1/3  and  the  impulse  sealing  relation  beeomes 


I  OC 


A,y2(R> 


A 


i  /2  /  r: 


a 


(1-12) 


where  H  is  the  total  energy  in  the  laser  pulse  and  is  the  pulse  length. 
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The  impulse  normalized  by  the  target  area  and  the  square  root  of  th  focal  spot  area  are 
presented  in  Fig.  1-21  as  a  function  of  the  average  power  density  as  suggested  by  Eq.  (1-11 ).  The 
solid  line  is  a  least-squares  fit  to  the  data  and  has  a  slope  of  0.27.  The  dashed  line  is  fit  to  the 
data  with  a  slope  of  1/3  as  predicted  by  Eq.  (1-11  ).  The  scatter  of  the  data  makes  it  impossible 
to  determine  the  precise  power  of  the  impulse  dependence  on  the  average  power  density,  but  a 
small  exponent  is  indicated  (i.e.,  l/5  to  l/3).  It  should  be  noted  that,  as  pointed  out  in  Ref.  1  3, 
the  power  density  is  only  known  to  within  about  a  factor  of  3.  Equation  (1-11  )  indicates  that  the 
normalized  impulse  is  relatively  insensitive  to  the  power  density  and  correlates  the  data  to 
within  about  a  factor  of  2. 


Fig.  1-21 .  Plot  of  normalized  impulse  as  a  function  of  average  power  density  for  data  obtained 
with  use  of  AERL  e-beam  COo  laser.  Solid  symbols  denote  AERL  data  from  Kef.  14.  Solid  line 
is  least-squares  fit  to  data,  1/Aj  =  2(q  )0. 27.  dashed  line  is  fit  to  data  with  slope  of  1/3. 

Increasing  the  total  energy  above  300  J  led  to  the  formation  of  multiple  aerosol-induced 
breakdowns  in  the  laser  beam.  An  open-shutter  photograph  of  this  effect  is  shown  in  Fig.  1-22. 
The  total  impulse  delivered  to  an  effectively  infinite  target  (40cm  )  was  measured  on  a  number 
of  runs  both  with  and  without  the  formation  of  aerosol- induced  breakdowns.  In  Fig.  1-2  3,  the 
measured  impulse  is  shown  as  a  function  of  total  energy  in  the  laser  pulse.  The  limited  data 
show  a  linear  increase  of  impulse  with  energy  in  agreement  with  Eq.  (1-12)  when  multiple  break¬ 
downs  were  not  observed,  and  an  apparent  saturation  of  delivered  impulse  when  multiple  break¬ 
downs  were  observed.  This  interpretation  of  the  data  indicates  that,  with  AERL  machine  pulse 
shape  (an  initial  large  spike),  the  increased  energy  above  300  to  400  J  is  being  absorbed  in  the 
breakdowns  and  is  decoupled  from  the  target*  j  p  jjOWCjer 

\u  iill>*niiiti\i'  i n 1 4-r|»ri*l ,it ion  11  f  these  il.it.i  is  pussilde  if  tin*  .issiiuiptiuns  uf  impulse  delixery  In  ,t  spherir.il  hlnst  n.ne 
f  ill  slum. mentis  ilr  jm- it  r.li  uf  enerpx  .it  .1  [mint}  .Mill  ,1  futile  l.trpel  .ire  lll.nlr.  11.12  'litis  ittnilel  prrdit  ts  llt.ll  the  impulse  is 
prnpurl  intnii  In  tiie  -ijii.if'e  runt  uf  tul.il  enerpx.  lite  impulse  il.it, ■  slinnti  in  |-'ip.  1-21  rnultl  he  fit  tu  uilliiti  experiment, tl  errur  In  <t 

fum  tieil  uf  the  re  runt  uf  the  enerpx .  This  inter |U  el  ill  inn  then  s.n  s  the  uernsnl-imlin  ell  lire, ltd,  mils  nhsnrli  ntiK  ,i  sin, ill  .it ml 

'f  ellerpl  , i tid  ilu  mil  .tffei  I  the  tiie.isured  i lit,  ttlse.  It  is  i  le.tr  th.it  mure  d.lt.i  .ire  net  ess.irx  tu  re. tils  settle  this  [mint. 
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Fig.  1-22.  Open-shutter  photograph  of  aerosol-induced  breakdowns  in  laser  beam 
which  is  focused  on  carbon  target  with  27-m  focal-length  mirror.  Total  energy  in 
this  shat  was  514  J.  ' 
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Fig.  1-23.  Plat  of  total  impulse  as  a  function  of  energy  in  laser  pulse.  Circles  denote 
observation  af  aerosol-induced  breakdowns  in  laser  beam.  Dashed  line  indicates  trend 
of  impulse  saturation  with  increasing  energy. 
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b.  Plasma  Properties  Associated  with  Irradiation  of  a  Surface  by  a  Pulsed  Laser 

A  number  of  experiments  were  performed  with  the  Avco  device  to  try  to  quantify  some  of  the 
properties  of  the  plasma  produced  in  the  laser-target  interaction.  The  results  of  some  earlier 
camera  data  (Ref.  11)  indicated  that  a  self-luminous  front  with  velocities  of  2  X  105  to  lO^cm/sec 
moved  back  in  the  direction  of  the  laser  beam.  To  supplement  this  information  and  provide  two- 
dimensional  data,  the  streak  camera  was  replaced  by  a  framing  camera  running  at  about  2  X  1  0^ 
frames/sec.  In  addition,  the  plasma  was  back-illuminated  by  a  40-psec-long  pulse  of  4881  radi¬ 
ation,  so  that  essentially  the  second  derivative  of  index  of  refraction  (shadowgraph)  rather  than 
the  self-lumination  was  viewed.  The  target  was  a  carbon  surface  and  the  incident  laser  energy 
was  approximately  200  J.  Some  typical  results  are  shown  in  Figs.  1-24  through  1-26, 

Figure  1-24  shows  the  growth  of  the  plasma  on  the  target.  The  growth  is  approximately  that 
of  a  spherical  wave,  with  some  asymmetry  resulting  because  the  impinging  laser  energy  tended 
to  arrive  in  two  separate  cones  rather  than  in  one.  Figure  1-25  shows  similar  data,  except 
that  in  this  case  a  fine  copper  wire  was  placed  in  front  of  the  target,  and  breakdown  initiated  on 
the  wire  shows  a  distinctly  spherical  expansion  with  an  expansion  velocity  of  1.8  x  lo5  cm/sec. 
Finally,  Fig.  1-26  shows  the  resulting  pattern  in  the  late  stages  of  expansion,  after  the  laser  pulse 
is  completed  and  the  initial  shock  wave  has  moved  away.  Notice  the  turbulent  state  of  the  medium. 
By  using  these  types  of  data  and  neglecting  the  asymmetry  due  to  the  laser  intensity  distribution, 
the  lateral  and  radial  growth  of  a  plasma  ball  was  measured  and  is  shown  in  Fig.  1-27, 

These  data  were  used  to  help  determine  the  ahsorptivity  of  the  plasma.  In  an  experimental 
arrangement  similar  to  that  described  in  Ref.  11  and  illustrated  in  Fig.  1-28,  the  transmission 
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Fig.  1-24.  Back-illuminated  framing  camera  data  showing  growth  of  a  plasma  on  a  target. 
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Fig.  1-25.  Bock-illuminoted  framing  data  showing  initiation  ond  growth  of  a  plasmo 
on  o  fine  copper  wire.  Interframe  time  is  0.5  psec. 
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Fig.  1-26.  Late  time  development  af  a  plasma  initiated  an  a  surface.  Nate  turbulent  nature 
of  plasma  medium.  These  data  are  typical  ~30  psec  after  initiation  af  10-pm  laser  pulse. 
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F'9'  'r27'  „Lateral  and  farword  growth  af  plasma  on  a  ccrban  surface,  fcnerqy  in  pulse 
was  about  200  J.  '  v 
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/l9-' i"D2^'jtiEXP!nment?  arran9erTlent  used  fa  measure  plasma  opacity.  An  ian-argan  laser 
(A  4881  A)  and  an  order  af  10-pm  high-power  laser  are  made  fa  traverse  same  path  thrauqh 
plasma  medium  a  distance  d  in  front  of  a  target.  M,  and  M2  are  bending  mirrors,  D,  and 

tendlm-t filte^09  defeC,arS'  °3  iS  °  phatadiode  dRtectar'  and  F,  is  a  chain  af  narraw- 
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history  of  10.6-gm  and  4881-A  radiation  through  the  plasma  was  monitored  a  distance  d  in 
front  of  the  target. 

The  10.6-pm  radiation  used  to  probe  the  plasma  was  obtained  by  taking  one  of  the  grating 
side  orders  and  transmitting  it  through  the  plasma  and  into  a  photon-drag  detector  D^.  At  the 
same  time,  the  unattenuated  laser  pulse  was  monitored  with  photon-drag  detector  D^.  A  1-W 
40-p.sec-long  pulsed  argon-ion  laser  was  triggered  approximately  lOpsec  before  the  C02  laser 
pulse  and  passed  through  the  same  path  as  the  deflected  side  order  of  the  C02  pulse.  The  argon 
laser  pulse  was  passed  through  a  series  of  4881-A  narrowband  filters  into  a  Si  photodiode  D3_ 


Fig.  1-29.  Normolized  10-pm  pulse  shapes 
o  different  distance  in  front  of  a  carban  tar¬ 
get.  Solid  curves  are  unattenuated  pulses, 
while  dotted  curves  represent  pulses  ofter 
being  transmitted  through  plasma  medium, 
(a),  (b),  ond  (c)  represent  data  taken  at  dis¬ 
tances  of  0,  3,  ond  10  mm  in  front  of  tar¬ 
get,  respectively. 


Hi)  I  .i  ta' 


Measurements  were  taken  at  three  distances:  0,  3,  and  10mm  in  front  of  the  target.  For 
the  10-p.m  data,  the  two  pulses  (attenuated  and  unattenuated)  were  normalized  to  peak  unit  inten¬ 
sity.  The  results  are  shown  in  Fig,  I-29(a-c).  Notice  that  the  C02  laser  pulse  has  a  rise  time 
of  400  to  500 nsec.  Absorption  begins  after  a  delay  which  is  consistent  with  velocities  in  the 
range  3  X  1 05  to  10^  cm/sec.  The  framing  camera  data  were  used  to  determine  the  optical  path 
length  <5(d,  t)  at  a  given  distance  d  in  front  of  the  target  and  at  a  given  time.  The  absorption 
coefficient 

'"'A, 

“  =-  ~6  — 


was  then  calculated  and  is  shown  in  Fig.  1-30.  The  simplifying  assumption  of  a  uniform  absorp¬ 
tion  across  the  plasma  was  used. 
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Fig.  1-30.  Time-dependent  absorption  coefficient  at  10.6  pm  at  three  distances 
in  front  af  target. 


Fig.  1-31 .  Argon-ion  laser  pulse  transmitted  through  absorbing  vapor  in  front  af  a  carbon  target, 
(a),  (b),  and  (c)  represent  distances  of  0,  3,  and  10mm  in  front  of  target,  respectively,  (d)  is  at 
0  mm  in  front  of  target,  but  one  4881 -A  filter  was  removed  from  front  af  detector  Dg  (Fig.  1-28) 
sa  that  self- luminosity  af  plasma  in  wavelength  regime  was  detected.  Notice  absorption  dip  that 
occurs  at  t  0,  before  self-luminosity  of  plasma  is  detectable. 
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Data  from  the  4881-A  laser  are  shown  in  Figs.  I-  3 1  (a - d ).  The  relative  transmission  at  d  0, 

3,  and  10mm  are  shown  in  Figs.  I-31(a-e),  respectively.  Figure  I- 31  (d )  is  taken  at  0mm,  but 
one  of  the  narrowband  filters  in  front  of  detector  D^  was  removed  so  that  the  plasma  self-emission 
is  shown.  Notice  that  there  is  a  narrow  absorption  spike  (<1  psec)  before  the  self-emission  of  the 
plasma  builds  up. 

By  analyzing  these  data  in  the  same  way  as  the  10-pm  data,  the  absorption  eoeffieient  at 
4881  A  is  shown  in  Fig.  1-32. 


Fig.  1-32.  Time-dependent  absorption  coefficient 
data  at  4881  A  at  three  distances  in  front  af  target. 
Because  af  jitter  in  firing  af  this  laser,  uncertainty 
in  initial  points  af  these  curves  is  ~1  psec. 
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If  one  makes  the  assumption  that,  for  the  10-pm  radiation,  the  absorption  is  due  to  inverse 
bremsstrahlung,  one  can  use  the  absorption  coefficient  to  determine  the  electron  density,15’^ 
using  the  formula 


rvfem'1)  =1.6  x  1  o'  37  A2  r;2  T  '  3^2  (1  -  e 

'e  e 


-hr/kT  kT 
e\  — § 
hr 


(M3) 


Here,  a  is  the  absorption  coefficient  in  cm'1,  A  is  the  wavelength  in  pm,  rjp  is  the  electron 
density  in  number/em3,  is  the  electron  temperature  in  eV,  and  v  is  the  laser  frequency. 

The  electron-density  histories  shown  in  Fig.  1-3  3  were  calculated  with  Kq.  (1-13),  the  absorption 
coefficients  given  in  Fig.  1-30,  and  the  assumption  of  a  1-eV  plasma.  Also  shown  in  Fig.  1-33 

o 

arc  electron  densities  obtained  in  the  same  manner  from  the  4881-A  absorption  data.  The  as- 
sumption  of  a  1-eV  plasma  temperature  is  a  reasonable  estimate  of  the  actual  temperature; 
thus,  the  calculated  electron  densities  are  of  the  right  order  of  magnitude  and  give  correct  rel¬ 
ative  values. 
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Fig.  1-33.  Electron  densities  calculated  from  dato  of  Figs.  1-30  ond  1-32  for  1-eV  temperoture 
ondassumption  thot  predominant  absorption  mechanism  is  inverse  bremsstrohlung.  This  ossumption 
is  plousible  for  10-pm  doto  and  gives  reasonable  estimote  for  electron  density  time  history  in 
10l  7  to  3  y  !0l  ^  cm '3.  The  4881  -  A  doto  suggest  some  photo-excited  absorbing  species  as  more 
proboble  condidates  for  vapor  opacity. 
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Those  data  indicate  that  the  absorption  at  4881  A  is  not  primarily  due  to  free  electrons.  In 
fact,  to  within  about  ~i  psec,  the  absorption  at  4881  A  begins  simultaneously  at  all  three  dis¬ 
tances  d  in  front  of  the  target.  Recall  the  ahsorption  dip  at  the  heginning  of  the  10-pm  laser 
pulse  in  Fig,  1- 31(d)  and  notice  that  the  ahsorption  takes  place  before  there  is  appreciable  self¬ 
emission  from  the  plasma,  i.e.,  the  plasma  is  still  relatively  cool  with  no  appreciable  electron 
density. 

The  simultaneous  nature  of  the  absorption  onset  suggests  the  production  of  a  pnoto- excited 
species  which  absorbs  strongly  at  4881  A  and  decays  with  a  lifetime  of  about  10  nsec. 

H.  Kleiman 
R.  W.  O'Neil 

3.  Explosive  Vaporization  of  Fog  Droplets  by  a  10.6-pm  Laser  Pulse 

Study  continues  of  the  exnlosive  vaporization  of  fog  droplets  by  a  10.6-pm  laser  pulse.  The 
early  stage  (up  to  about  2 psec)  of  the  explosion  process  was  described  in  previous  Optics  Re¬ 
search  Reports.*  We  studied  the  later  stage  of  the  process  by  observing  the  vapor  "blob"  pro¬ 
duced  by  the  evaporation  of  single  droplets.  These  observations  were  made  with  a  high-speed 
schlieren  photography  system  that  was  described  in  the  last  Optics  Research  Report  [(1971:2), 
DDC  AD-90121  3);  early  results  on  the  vaporization  of  oversize  droplets  were  also  presented  in 
that  report. 

In  a  ground  fog  or  a  cloud,  the  water  droplets  usually  range  in  size  from  <1-  to  ~25-pm  ra¬ 
dius.  In  the  present  study,  droplets  ranging  from  5-  to  25-pm  radius  have  been  irradiated.  In 
addition,  some  oversize  droplets  (25-  to  70-pm  radius)  have  been  studied  because  their  behavior 
does  help  in  understanding  the  processes  involved.  The  laser  energy  at  the  position  of  the  water 
droplet  was  varied  between  5  and  15  J/cmZ  in  these  experiments.  The  peak  power  density  was 
typically  10  to  30  MW/cmZ. 

Electrical  breakdown  did  not  occur  under  these  conditions.  This  was  determined  bv  means 
of  time-exposed  photomicrographs  of  the  exploding  droplet;  these  photos  were  taken  without  any 
back-illumination.  No  light  was  emitted  during  the  explosive  vaporization  of  a  water  droplet. 

If,  however,  instead  of  a  water  droplet,  a  25-pm-radius  p  ass  fiber  is  placed  at  the  same  posi¬ 
tion  and  irradiated  with  an  identical  pulse,  electrical  breat%  'own  does  occur  and  a  brilliant  flash 
of  light  is  produced. 

a.  Results  and  Discussion 

The  schlieren  pictures  show  an  expanding  sphere  of  vapor  surrounded  by  a  rapidly  expanding 
shock  wave.  Actually,  the  observed  "vapor  sphere"  is  a  sphere  of  air  that  has  been  heated  hy 
the  vapor.  The  sphere  of  heated  air  has  an  index  of  refraction  different  from  that  of  the  ambient 
medium,  and  this  difference  is  sensed  by  the  schlieren  svstem. 

Figure  1-34  shows  schlieren  photographs  of  four  exploding  droplets  at  times  from  0.7  to 
37  psec  after  the  10.6-pm  laser  pulse.  The  radius  of  each  water  droplet  is  indicated  in  the  pic¬ 
ture.  At  0.7-psec,  the  vapor  sphere  cannot  be  distinguished  within  the  shock  wave,  but  the 
second  photograph  (4.9  psec)  clearly  shows  the  vapor  sphere  at  the  center  of  the  shock-wave 

(►pill  -  Ilf- iMirli  lli.|iurl-.  1  irri  oln  I ..ilmr.it nr a  .  M.l.T.  (1970:.'!).  I)I)C  AD-IWJO 1 7,  Hn(|  (1971:11.  DDC  \D-HHKH2k-l .. 
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Fig.  1-34.  Schlieren  photographs  of  four  exploding  woter  droplets.  Rodius  of  eoch  droplet 
is  shown  ot  left,  ond  elopsed  time  since  10.6-pm  loser  pulse  is  shown  ot  right.  Loser  pulse 
direction:  left  to  right. 

sphere.  In  the  third  picture  (16.2  psec),  the  shock  wave  can  be  seen  a'  he  top,  after  it  has  re¬ 
flected  from  the  horizontal  surface  (which  happens  to  be  the  bottom  of  the  droplet  generator). 
Finally  (at  37  psec),  the  shock  wave  is  beyond  the  FOV  of  the  system.  This  sequence  of  pictures 
also  shows  the  expanding  volume  of  the  vapor  sphere.  In  these  photographs,  the  direction  of  the 
10.6-pm  laser  pulse  is  from  left  to  right;  the  schlieren  knife-edge  is  oriented  horizontally.  (The 
pictures  also  show  some  extraneous  features  that  are  common  to  all  the  photographs;  these  are 
due  to  aberrations  in  the  optical  system  and  to  imperfections  in  the  lenses.  They  can  be  readily 
distinguished  and  disregarded.) 

In  the  study  of  the  early  stage  of  the  water-droplet  explosion,  we  showed  that  the  liquid-vapor 
mixture  was  expanding  at  a  radial  velocity  of  0.13mm/psec  in  the  time  interval  0  to  2.  3  psec. 

Our  schlieren  photographs  show  that  the  vapor  sphere  expands  with  a  radial  velocity  of  0.015  mm/ 
psec  in  the  time  interval  5  to  16  psec,  and  the  velocity  has  decreased  to  a  value  of  O.OlOmm/psec 
in  the  interval  16  to  37  usee.  Kven  in  this  interval,  the  radial  expansion  rate  is  ~300  times 
greater  than  the  diffusion  velocity,  indicating  that  the  expansion  process  is  still  dominated  by 
the  momentum  of  the  initial  explosion. 

In  Fig.  1-35,  the  photographs  show  the  effect  of  droplet  size  on  the  magnitude  of  the  ohserv- 
ahle  effects.  The  small  droplets  (11-  and  13-pm  radius)  at  the  left  produced  a  faint  but  clearlv 
observable  shock  wave  with  a  faintly  distinguishable  vapor  sphere  at  the  center.  The  pictures 
at  the  right  show  the  much  stronger  effects  produced  by  the  explosive  vaporization  of  25-pm- 
radius  droplets,  where  the  mass  of  water  is  ~8  times  as  great.  In  each  of  these  shots,  the  energy 
density  of  the  laser  pulse  was  about  10  j/cm^,  and  the  pulse  shape  was  about  the  same. 
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Fig.  1-35.  Schlieren  photographs  of  four  exploding  woter  droplets,  showing  effect  of  droplet 
size.  Loser  pulse  direction:  left  to  right. 


Under  certain  conditions,  the  hot  vapor  assumed  a  nonspherical  shape  -  usually  when  the 
water  droplet  was  very  large  (> 35-pm radius) or  when  the  energy  density  was  very  low  (~5  j/cm2). 
In  these  cases,  the  vapor  volume  was  shaped  like  a  dumbbell  with  its  axis  oriented  approximately 
along  the  direction  of  the  laser  pulse.  Figure  1-36  shows  a  typical  example,  an  ~35-pm-radius 
droplet  at  6.7  psec  after  the  laser  pulse.  The  left  sphere  of  the  dumbbell  is  at  the  center  of  the 
shock  wave,  and  this  is  the  site  of  the  original  explosion.  If  the  absorbed  energy  is  not  sufficient 
to  cause  immediate  (within  ~2  psec)  vaporisation  of  the  droplet,  and  if  the  droplet  is  large  enough 
to  undergo  surface  heating  rather  than  volume  heating,  then  the  remainder  of  the  evaporating 


Fig.  1-36.  Schlieren  photograph  of  lorge 
(>  35-pm  radius)  exploding  water  droplet 
6.7  psec  after  10.6-pm  laser  pulse.  Loser 
pulse  direction:  left  to  right. 
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droplet  is  propelled  by  the  rocket  action  of  its  vaporizing  surface.  The  motion  is  always  oppo¬ 
site  to  the  direction  that  the  blow-off  vapor  is  moving,  as  required  to  conserve  momentum. 

A  picture  similar  to  Fig.  1-36  was  obtained  with  a  26-pm-radius  droplet  irradiated  with  a 
laser  pulse  of  4  J/cm2.  But  when  another  droplet  of  the  same  size  was  irradiated  with  a  9-J/cm 
pulse,  the  usual  spherical  vapor  volume  resulted.  This  behavior  can  be  explained  by  the  curves 
in  Fig.  1-37,  which  show  the  energy  absorbed  by  a  water  droplet  for  pulses  of  various  energy 
densities.  Actually,  this  is  the  energy  that  would  he  absorbed  by  a  water  droplet  that  remained 
into et  during  th'-  laser  pulse.  These  curves  were  drawn  using  the  absorption  cross  sections  in 
the  tables  of  Zel'manovich  and  Shifrin.17  In  the  case  of  the  9-J/crr/  pulse,  the  energy  that 


Fig.  1-37.  Energy  absorbed  by  woter  droplet  for  incident  10.6-pm  laser  pulses  of  indicoted 
energy  density.  Energy  of  voparizatian,  E^  ,  's  °^so  shown  for  comparison. 

could  be  absorbed  by  the  water  droplet  is  substantially  greater  than  the  energy  required  for 
vaporization.  In  the  case  of  the  4- J/cm2  pulse,  the  energy  absorbed  by  the  droplet  is  much  less 
than  the  energy  of  vaporization.  Our  early  studies  using  the  tri-scopic  camera  system  showed 
that  the  droplets  of  radius  >  15  pm  undergo  surface  heating.  Such  droplets  when  partially  vapor¬ 
ized  can  be  propelled  along  the  general  direction  of  the  laser  pulse;  this  is  the  case  for  the  drop¬ 
let  in  Fig,  1-  16,  which  could  absorb  barely  enough  energy  for  vaporization.  In  general,  the 
asymmetric  heating  (surface  heating)  of  optically  thick  droplets  will  result  in  an  asymmetric 
vapor  volume  unless  the  absorbed  energy  is  much  greater  than  that  required  for  vaporizing  the 
droplet.  The  asymmetric  volume  will  be  shaped  like  a  dumbbell  or  like  an  incipient  dumbbell, 
with  the  second  sphere  undeveloped  or  partially  developed. 
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b.  Energy  of  the  Explosion 

Figure  1-38  shows  the  radius  of  the  shock  wave  as  a  function  of  time  for  manv  different 
shots,  covering  three  ranges  of  droplet  size.  Although  the  points  fall  close  to  the  line  of  slope 
corresponding  to  a  shock-wave  velocity  of  Mach  1  (in  air),  it  will  be  noted  that  the  zero-time 
intercept  is  0.5  mm,  which  is  20  times  greater  than  the  initial  droplet  radius.  In  the  first  couple 
of  microseconds,  the  shock-wave  velocity  is  much  greater  than  Mach  1. 
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Fig.  1-38.  Radius  of  shock  wove  as  function  of  time. 

The  radius  of  the  shock  wave  at  a  given  time  can  be  used  to  calculate  the  explosive  energy, 
as  discussed  in  the  last  Optics  lie-search  lieport  [(1971:2),  1)1)0  AD-901213,'.  We  have  calculated 
the  explosive  yield  for  5  shots  involving  water  droplets  of  25-pm  radius.  The  yield  may  be  com¬ 
pared  with  the  difference  between  the  absorbed  energy,  as  determined  from  Fig.  1-37,  anti  the 
energy  needed  to  evaporate  the  droplets.  The  agreement  is  within  about  25  percent.  This  en¬ 
ergy  balance  cannot  be  used  for  large  (  >40-nm-radius)  droplets  because  they  do  not  evaporate 
completely. 

P.  Kafalns 
J.  Ilerrmann 
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Fig.  1-39.  Schematic  diagram  of  500-J  EDL. 
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C.  DEVICES 

1.  Electron  Beam  (e-Beam) 

a.  500-J,  e-Beam  Excited  Laser 

Figure  1-39  is  a  schematic  diagram  of  the  500-J  EDL  now  under  construction.*  All  major 
components  are  on  hand  and  assembly  is  nearly  complete.  In  ke  ping  with  its  planned  use  as 
a  source  of  high-energy  10,6-p.m  pulses  of  good  optical  quality,  the  operating  sequence  has  been 
automated  as  much  as  possible  to  keep  the  output  reproducible  and  reliable. 

The  electron  gun  is  a  triode  and  can  be  grid  controlled  from  cutoff  to  current  densities  in 
2 

excess  of  lOOmA/cm  through  a  15  X  104cm,  1-mil  titanium  foil  window.  It  is  capable  of  pulses 
of  from  1-  to  100-psec  duration  with  rise  times  of  0.1  psec. 

The  optical  cavity  has  an  effective  length  of  2.5  m  with  a  maximum  cross  section  of  15  X  20cm. 
The  active  volume  is  22.5  liters. 

To  avoid  the  need  for  large,  high-power  windows,  a  fast  shutter  valve  is  opened  to  pass  the 
outcoupled  optical  beam. 

The  optical  cavity  will  be  evacuated  after  each  shot  and  then  refilled  with  the  laser  gas  mix¬ 
ture.  A  gas  blender  will  allow  the  mixture  to  be  easily  optimized  for  the  pulse  length  and  power 
chosen  by  an  experimenter. 

First  operation  in  a  laser  is  expected  early  in  August,  followed  by  about  60  days  of  tests 
and  diagnostics.  We  anticipate  that  the  device  should  be  operational  in  October. 

A.  J.  Morency 

R.  J.  Carbone 

b.  e-Beam  C02  Laser  Kinetics 

1 8 

The  C02-N2~He  rate  equations  of  Douglas- Hamilton  were  programmed  using  the  Runge- 
Kutta  method.  In  Douglas-Hamilton's  approach,  the  v j,  and  modes  of  CO.,  and  the  vibra¬ 

tional  levels  of  N2  are  assumed  to  equilibrate  internally  so  rapidly  that  each  of  these  modes 
can  be  described  by  a  vibrational  temperature.  The  four  equations  describing  these  temperatures, 

plus  a  fifth  equation  describing  the  gas  temperature,  compose  the  coupled  set  of  differential  equa- 

19 

tions.  With  the  electron  excitation  coefficients  calculated  by  Nighan  and  the  relaxation  rates 
reported  by  Taylor  and  Bitterman,^0  the  computed  gains,  pulse  lengths,  and  pulse  energies  were 
in  reasonable  agreement  with  the  experimental  results  on  the  Avco  e-beam  laser. 

In  this  analysis,  the  laser  flux  is  determined  by  the  constraint  that,  during  the  pulse,  the 
laser  gain  is  equal  to  the  cavity  losses.  Accordingly,  these  equations  do  not  correctly  describe 
the  initial  portion  of  the  pulse  where  the  gain  still  greatly  exceeds  the  losses,  and  the  equations 
do  not  predict  the  experimentally  observed  initial  spike.  To  remedy  this,  a  sixth  equation  de¬ 
scribing  the  growth  of  the  stimulated  emission  was  included.  This  indeed  provided  the  initial 
spike,  but  its  computed  50-nsec  width  is  considerably  smaller  than  the  observed  300  to  500 nsec. 
This  discrepancy  may  be  due  to  the  fact  that  on  the  time  scale  of  the  initial  spike,  the  assump¬ 
tion  of  vibrational  temperatures  may  be  invalid.  We  are  presently  investigating  this  possibility. 

S.  Marcus 
D.  Zanni 

•  >r,-  OplKs  U,..s<-aril,  lirpnrl  (1971:2),  p.  23,  l)l)C  AD-901213. 
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Fig.  1-40.  Geometry  of  500-J  pulse  unstoble  confocol  resonotor  covity  for  optimum  mode 
quol  ity-energy  product  (outcoupling  S  =  0.85). 


GEOMETRIC  OUTPUT  COUPLING,  8 

Fig.  1-41  .  Output  intensity  normalized  to  soturotion  intensity  os  function  of  geometric  outcoupling 
(5  |  _  (i/fn  2)]  for  uniform  mode  filling,  with  mirror  reflectonce  =  0.99,  covity  length  2.5  m, 

octive  length  1.0  m,  ond  smoll  signal  goin  0. 032  cm"^ . 
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c.  Cavity  Design  of  a  Pulsed  500-J  CO.,  Laser 

As  its  first  priority,  the  EDL  effort  has  a  pulsed  CO^  laser  device  that  is  designed  to  out- 
couple  500  J  of  10.6-nm  radiation.’’''  The  volume  of  the  laser  medium  will  be  22.5  liters  con¬ 
tained  within  a  length  of  1  m.  A  2.5-m  unstable,  confocal,  positive  branch  resonator  (shown  in 
Fig.  1-40)  was  chosen  to  effectively  utilize  this  large  volume/short-length  gain  medium. 

The  cavity  design  was  chosen  so  that  optical  energy  will  not  focus  within  the  device,  and  the 
forward  and  backward  energy  flow  will  very  nearly  fill  the  gain  medium.  The  optimized  reso¬ 
nator  parameters  were  chosen  as  a  tradeoff  between  total  power  output  reduction  while  maintain¬ 
ing  a  single  mode  as  the  output  coupling  was  increased.  Figure  1-41  illustrates  the  effect  of 

-1  21 

power  output  vs  coupling  at  a  gain  of  0.032  cm  and  99-percent  reflecting  mirrors.  The  Rigrod 

approach  was  used  to  determine  the  power  extraction  as  a  function  of  predicted  gain  and  optimal 

energy  storage.  A  reasonable  tradeoff  between  mode  quality  and  optimum  energy  output  is  at  an 

output  calculated  to  be  60  percent.  Using  a  geometrical  approach,  this  outcoupling  will  utilize 

at  least  52  percent  of  the  active  volume  for  an  inscribed  (within  a  15  x  15  cm  square)  circular 

mode,  65  percent  for  a  15  x  15  cm  square  mode,  and  87  percent  for  a  fitted  rectangular  mode  of 

15  x  20cm.  Several  cavity  outcouplings  will  be  tested  to  experimentally  verify  the  power  output 

predictions,  and  to  compare  the  experimental  results  with  inode  calculations  that  have  been  re- 
22-25 

ported  and  with  calculations  presently  being  done  at  Lincoln  Laboratory. 

— Thc^jf>H?if-to[rL.'ffPcrT?f',fudge  beem  quality  outcoupled  from  the  laser  is  the  energy  deposited 
in  the  Airy  disks  compared  with  81  percent  of  the  energy  deposited  in  the  first  Airy  disk  for  a 
uniformly  illuminated  aperture.  For  the  60-percent  coupling  case,  83  percent  of  the  output  en¬ 
ergy  will  be  deposited  within  the  first  three  Airy  disks,  or  3X  diffraction-limited  by  this  criterion. 

Optimization  of  mode  quality  is  achieved  by  maximizing  the  product  of  the  energy  output  and 
the  energy  fraction  in  the  first  Airy  disk,  as  shown  in  Fig.  1-42.  This  optimum  outcoupling  is  at 
85  percent,  which  corresponds  to  a  reduction  of  about  15  percent  in  energy  extracted  per  unit 
volume  (relative  to  60-percent  outcoupling),  approximately  82-percent  volume  filling  factor, 
and  about  45  percent  of  the  net  energy  within  the  first  Airy  disk. 

The  absolute  power  outcoupled  for  the  cavity  design  has  been  calculated  from  the  rate  equa¬ 
tions  and  energy  balance  within  the  laser  medium.  The  approach  used  followed  that  of  Douglas- 
1 8 

Hamilton.  Predictions  of  medium  heating,  pulse  length,  and  gain  are  reliable,  but  calculated 

2  2 

optical  energy  storage/cross-sectional  area  at  7.4  j/cm  with  2mA/cm  e-beam  pumping  may 
be  too  high.  This  indicates  an  extraction  of  74  J/liter  for  aim  device.  Table  1-1  gives  typical 
computed  values  of  pulse  length,  gain,  and  output  energy  for  different  e-beam  current  intensities. 

R.  J.  Carbone 
H.  Granek 


d.  Mode  Calculations  for  the  Unstable  Resonator 

The  Fredholm  integral  equation  for  the  asymmetric  unstable  resonator  cavity  is  ir:ed  for 

the  mode  field  distributions  without  a  gain  medium  in  the  mirror  planes.  The  field  distributions 
22 

are  given  by 


'  Optics  Hcsc.iri  li  llcpnrl  (!l171:2),  11.  git.  !)!)(!  \I)-U0)2M. 
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Fig.  1-42.  Amount  of  output  intensity  (normolized  to  soturotion  intensity)  concentrated  within 
first  Airy  dork  ring  of  for-field  pattern  (incident  plone  wove  opproximotion  onto  outcoupling 
annulus)  for  conditions  given  in  Fig.  1-41. 


TABLE  1-1 

TYPICAL  COMPUTED  LASER  PARAMETERS  FOR  GIVEN  e-BEAM 
CURRENT  INTENSITIES 

Current  Density 
(mA/cm2) 

Peok  Zero-Flux  Gain 
(cm-1) 

Pulse  Length 
(psec) 

Energy  Out 
(J/cm2) 

2 

0.0282 

43.8 

wm 

10 

0. 0326 

25.7 

20 

0.0339 

20.1 

1 

30 

0.0345 

17.4 

m 
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VlVl(x)  =  ji+1(k/L)  (  yJ/( kxy/L)  exp[-j(k/2L)  (gtx2  +  g2  y2)]  <p2(y)  dy  (1-14) 

yZ9Ziy)  =  Ji+1<k/L>  1  (kyz/L)  exp{-j(k/2L)  (g2y2  +  g^2)]  ^(z)  dz  .  (I-15) 

In  this  expression,  g.  =  1  -  (L/R.),  where  R.  is  the  radius  of  curvature,  L  isthe  mirror  separa¬ 
tion,  y's  are  eigenvalues  to  be  determined,  and  2a.  is  the  mirror  diameter. 

One  can  substitute  Eq.  (1-15)  into  Eq.  (1-14)  to  determine  an  eigenvalue  equation  for  the  field 
at  the  output  mirror.  Use  of  Bessel  function  identities2^  permits  elimination  of  the  integral  of 
the  product  Bessel  functions  and  reduction  to  a  form 


viv2”i(x) = fo  1  e*p[<-  <2glg2  -D(4|jL)j 
XJf  (z^L  xz)  z«V*>dz  ■ 

Manipulation  of  this  expression  and  expansion  of  the  unknown  functions  in  the  form23 

*W.«p[-i^(2gl62-.]x2]  =  £  «mU''(x2| 


m=0 


2 

where  Cj  =  2nai  /L  allows  simplification  to  an  infinite  matrix  of  the  form27 


(1-16) 


(1-17) 


(VlV 


(l  +t),' 

— n 


exp{-U-  [2j(2glg2  -  1)  +1]}  (U-)/+t 


x  Z  a. 
m=0 


L(/)  (V) 
mm  '  ' 


(1-18) 


Calculations  of  a  truncated  form  of  this  matrix  for  the  empty  cavity  modes  of  the  unstable  res¬ 
onator  are  under  way.  Comparison  of  these  results  with  propagation  calculations  and  experi¬ 
mental  data  is  anticipated. 

H.  Granek 


e.  Longitudinal  e- Beam  Excited  Laser 

The  3-nsec,  600-keV,  e-beam  source  described  in  the  previous  Optics  Research  Report 
was  used  to  ionize  a  CO^-He  medium  from  which  laser  pulses  lasting  up  to  120psec  were  ob¬ 
tained.  In  the  configuration  being  used,  the  e-beam  is  nearly  colinear  with  the  laser  axis  and 
normal  to  tte  sustainer  electric  field.  An  article  describing  the  operation  of  this  laser  has  been 
published.  In  the  meantime,  the  maximum  output  energy  was  increased  to  22,4  J  over  the 
11.5  J  reported  in  that  paper.  This  was  accomplished  by  improving  the  design  of  the  electrodes 
and  by  increasing  the  active  laser  volume.  We  are  presently  redesigning  the  laser  cavity  to 
permit  passive  Q-switching. 

S.  Marcus 
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f.  Pulse  Shaping  a  CC^  I.aser  Beam 

For  a  number  of  applications  it  is  desirable  to  be  able  to  shape  in  time  high-energy  laser 

pulses  of  10.6-pm  radiation.  One  approach  is  to  try  to  build  pulse-shaping  capabilities  into  the 

laser  itself,  for  instance,  hy  modulating  the  grid  of  an  e-beam  pumped  laser.  However,  for 

other  classes  of  pulsed  machines  it  is  more  practical  to  try  to  develop  external  pulse- shaping 

capability.  Here,  we  report  some  successful  preliminary  attempts  to  generate  simple  waveforms 

using  saturable  absorbers  external  to  the  laser  cavity.  A  similar  technique  has  recently  been 

29 

reported  for  use  in  the  visible  range. 

The  basic  idea  can  be  understood  by  considering  the  expression  for  absorption  in  a  two-level 
system: 

ri  +  f  exp  {-  [1  +  (i/is)]t/T}-| 


*■  s  J 

Here,  <v  is  the  unsaturated  absorption  coefficient,  I  =  hr/2crTis  the  saturation  intensity,  a  is 
the  cross  section  for  molecular  absorption,  i>  is  the  laser  frequency,  and  r  is  the  relaxation 
time  of  the  upper  level.  In  all  useful  absorbers,  the  initial  and  final  levels  interact  strongly  with 
closely  lying  rotational  and  vibrational  levels  and,  consequently,  the  two-level  picture  is  not  a 
valid  one.  However,  it  can  be  shown  that,  when  the  relaxation  time  among  these  levels  is  short 
compared  with  t  and  the  laser  pulse  length,  the  above  relation  still  holds  except  that  I  =  hv/2ar6, 
where  the  factor  l/d  is  the  effective  number  of  interacting  levels. 

This  gives  the  medium  saturation  time  constant 


T 

S 


Suppose  we  require  the  pulse  length  =  Tg.  This  would  allow  the  medium  to  saturate 
gradually  and  generate  a  ramp  function,  which  would  require  an  intensity 


/  T  , 

\  _  hr  /  t  ,\ 

(-  ~i 

)  Zard  It  ) 

'  P 

/  V  p  / 

If  the  peak  of  the  ramp  is  to  be  substantially  unattenuated,  then  we  must  have  I  »  I  or, 
equivalently,  r  »  r  .  In  this  case,  a  square  pulse  of  intensity 

I  =  ,  X.  =  -■>— 

ST  Zut  6 
p  p 

should  generate  a  ramp  function.  Figure  1-4  3  shows  the  experimental  arrangement  in  which  we 
used  a  CO^  laser  which  generated  a  pulse  with  uniform  power  output  of  s;4  X  to'’ W  for  30psec 
(see  lief.  28).  This  pulse  was  focused  through  a  20-cm  absorption  cell  with  Brewster  angle  win¬ 
dows.  A  KC1  flat  diverted  about  5  percent  of  the  energy  to  a  BaTiO^  calorimeter  so  that  the 

pulse  energy  could  be  monitored.  The  pulse  at  the  cell  exit  was  monitored  with  a  photon-drag 

2 

Ge  detector.  The  beam  was  approximately  1  cm  in  area  at  the  entrance  and  exit  of  the  cell,  and 
perhaps  l/4  that  area  at  the  center  of  the  cell. 
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|u-5-40H| 


Fig.  1-43.  Experimental  arrangement  for  observing  pulse  shaping. 

Figure  I-44(a-d)  shows  a  series  of  waveforms  generated  at  different  pressures  of  SF^.  In 
the  region  where  a  is  a  function  of  pressure,  the  pressure  and  intensity  can  be  adjusted  to  main¬ 
tain  the  product  Iff  constant.  This  allows  generation  of  the  same  waveform  at  varying  intensity 
levels. 

Figure  I- 45 (a -  d )  shows  the  results  of  using  Freon  12  as  an  absorbing  gas.  The  relaxation 
time  here  is  short  compared  with  the  pulse  length,  so  that  the  intense  leading  edge  of  the  pulse 
goes  through  ur.attenuated  while  the  rear  of  the  pulse  is  absorbed.  This  is  the  more  familiar 
case  of  pulse  shaping. 

We  have  demonstrated  the  ability  to  produce  slow  or  fast  opening  attenuators  for  uniform 
laser  pulses,  and  to  use  these  attenuators  to  produce  ramp  and  spikes  of  desired  temporal 

characteristics.  ,r  ,,,  . 

II.  Kleinian 

S.  Marcus 


g.  Feasibility  Study  for  Amplification  of  Short  Pulses 
and  Frequency- Modulated  Pulses 

Here,  we  report  the  results  of  a  feasibility  study  on  the  design  of  a  very  broad-band,  high- 
power  laser  amplification  system.  There  is  current  research  interest  in  an  amplifier  at  10.6  pm 
that  could  handle  either  a  repeated  burst  sequence  of  nanosecond  pulses  or  a  repeated  sequence 
of  pulses  less  than  a  microsecond  in  length  with  FM  up  to  1.5  0117..  Kncrgy  outputs  in  excess  of 
3  J  per  pulse  and  average  power  of  20  kW  are  required,  with  total  coherence  time  of  the  pulse 
sequence  extending  from  circa  30psec  to  5  msec.  To  handle  the  information  bandwidth  of  the 
pulses,  the  amplifier  would  have  to  operate  at  pressures  of  at  least  150  to  200  torr  for  500-MIIz 
signals,  and  approaching  an  atmosphere  of  pressure  for  full  bandwidth  capacity.30"  32  There  are 
four  devices  in  the  "high" -pressure  regime  Hiat  could  possibly  satisfy  the  requirements  of  min¬ 
imum  optical  distortion  (<k/40):  longitudinal  discharge,  cross-flow  laser33;  waveguide  laser34- 
preionized  devices;  and  independently  ionized  devices.  The  first  type  has  been  limited  to  pres¬ 
sures  below  100  torr,  thus  does  not  exhibit  sufficient  bandwidth  to  be  applicable.  The  second  is 
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(c)  (d) 


Fig.  1-44.  Pulse  shaping  using  SF^  in  absorption  cell.  Pressure  of  SF$  is  (a)  0,  (b)  10, 
(c)  20,  and  (d)  35  torr,  respectively. 


■H 


Section  I 


[  -5-3S871 


(C)  (d) 

Fig.  1-45.  Pulse  shaping  using  Freon  12  in  absorption  cell.  Pressure  of  Frean-12  is  (a)  0, 
(b)  40,  (c)  100,  and  (d)  200  tarr,  respectively. 
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too  small  fon  appreciable  energy  extraction  (although  potentially  useful  in  preamplification).  In 

35  36 

the  third  category,  three  electrode  devices  ’  are  a  possibility,  but  the  discharge  optically 

distorts  the  medium  (especially  in  CO  -rich  mixtures)  and  the  gains  are  appreciably  lower  than 

37  38 

in  the  last  category.  However,  the  e-beam  preionized  laser  ’  (see  also  Sec.  O-l-a)  is  a  de¬ 
vice  that  already  approaches  the  capabilities  necessary  for  the  application.  Likewise,  the  UV 
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photoionized  laser  presently  being  researched  (see  also  Sec.  C-2  a)  offers  similar  promising 
possibilities  for  this  application. 

lloth  in  the  pulse  burst  and  FM  pulse  modes,  the  problems  of  medium  homogeneity  and 
preservation  of  gain  (and  pulse  coherence)  during  variable  time  intervals  are  important,  arid 
shall  be  treated  subsequently.  Likewise,  temporal  pulse  distortion  in  both  cases  is  important. 
However,  there  are  differences  that  must  be  treated  separately,  the  ability  of  the  molecular 
medium  to  be  reasonably  utilized  in  times  short  compared  with  its  relaxation  rates,  spurious 
FM  effects  due  to  medium  dispersion,  saturation,  and  pulse  sharpening. 

The  following  relaxation  time  calculations^11  were  done  fora  IleiN^iCO^  mixture  of  3:2:1  at 

a  total  pressure  of  1  atm.  The  results  are  easily  scalcablc  in  pressure.  The  inter-rotational 

4 1 

relaxation  rate  is  ~0. 12  nsec-atm  ;  therefore,  sufficiently  strong  nanosecond  pulses  will  dump 
all  the  rotational  levels  of  the  001  vibrational  state,  while  such  a  pulse  would  be  too  fast  to  dunip 
any  of  the  higher  vibrational  00r^  states  of  CO.,  [r00  >  3  nsec-atm  (CO.,)]  or  allow  replenish¬ 

ment  by  the  nitrogen  pump  (rv  ~  120nsec-atm).  Since  these  two  mechanisms  contribute  to  the 

1  2 

efficiency  of  energy  extraction  in  CO  ,  the  saturation  intensity  would  be  reduced  by  at  least  an 
order  of  magnitude  for  individual  pulses.  However,  the  higher  vibrational  states  woidd  partially 
repump  before  the  next  pulse  comes  through,  and  the  nitrogen  would  weakly  repump  during  a 
sequence  of  such  pulses. 

A  single  arbitrary  pulse  (in  time)  passing  through  a  high  gain  amplifying  laser  medium  will 
distort  in  shape,  because  its  leading  edge  always  secs  an  unsaturated  gain  and  will  grow  expo¬ 
nentially  (until  other  loss  mechanisms  c.ome  into  play),  wdiile  the  back  sections  of  the  pulse  will 
see  a  successively  more  saturated  gain  medium.  This  results  in  pulse  sharpening  at  the  leading 
edge  (and  a  corresponding  frequency  domain  distortion  that  is  important  in  the  FM  mode).  For 
the  nanosecond  pulses,  the  CO?  laser  acts  as  a  two-level  (total  vibrational  state  populations) 

system  and  the  interaction  of  the  population  inversion  and  optical  field  growth  can  be  shown  to 
•12 

be  described  by 

Po(x)  exp [-Na  f*  p0(x)  dx] _ 

i;(x,  t)  -  -  “TT'v/c  —  ~~ 

exp  | Nffc  f  ‘ /L  V>0(t' )  dt']  +  exp|-\a  |n  po(x')  dx1]  1 


and 


<ro(t  -  x/c ) 

r(X,  t)  -  -  f|_y  /r'\ 

4  —  { 1  —  exp[~Na  7)o(x')  dx']}  exp  | -Mere  J  '  ' 


e>f)(t' )  dt'  ] 


where 


^  n .  +  n 
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and 


<!>  number  of  photons/volume 

m  =  particle  density  in  respective  vil)rational  states 
0  =  resonance  absorption  cross  section/active  molecule. 

For  a  uniform  axial  gain  profile  r)(x,  0)  =  ?j  and  a  specified  pulse  form  >p(0,  t),  the  output  pulse 
can  then  be  integrably  determined.  The  distortion  of  a  I/srentzian  pulse  of  shape 

...  2I)T  /  i  \ 

^in  1  '  NtrC  (t2  +  T2j 

has  an  output  shape  of 


out 


_ (21TT/Nitc)  |l/(tZ  FT2)] _ _ 

1  -  f  1  -  exp  exp  {-Da  |t  +  (2/ir)  tan  ^(t/T)]} 


as  shown  in  Fig.  1-46. 


Fig.  1-46.  Photon  pulse  shape  of  output 
of  amplifier  normolized  to  input  for  a 
Larentzian-shoped  pulse.  I  Reprinted  with 
permission  from  E.L.  Steele,  Opticol  Lo¬ 
sers  in  Electronics  (Wiley,  New  York, 
1968),  p.  158.] 


A  square  pulse  going  through  the  medium  experiences  a  much  more  severe  sharpening  at 
t he  front  edge: 


24- 


nT_1  (t>  -  u.jU-  to)j 


(2<|,„/N)  -  U _ , ( t  —  to)J 


<'out^  1  -  (i  —  exp  [-  tr  1/ 1 . 1 }  exp(-20cl'otj 


as  exhibited  in  h’ig.  1-47. 

Independen'  of  the  particular  waveform,  we  can  estimate  the  total  energy  that  a  p  ise  of  a 
given  input  energy  will  extract  from  the  medium,  with  a  knowledge  of  the  small  signal  gain  and 
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Fig.  1-47 .  Ratio  of  output  photon  density 
H|_(t)  to  input  photon  density  qQ  vs  time  in 
units  of  pulse  width  t  for  on  initiolly 
squore  pulse  troversing  on  omplifying  me¬ 
dium  of  length  L=  10  cm,  initial  differ¬ 
ence  of  excited  ond  ground  stote  densities 
^o  ~  8  ^  10'8  cm-2,  ond  resononce  cross 
section  a  =  2.5  X  10-20cm2.  Totol  num¬ 
ber  of  photons  per  unit  areo  in  pulse  is 
q  =  4  X  10  cm-2  for  curve  o,  ond  q  = 
4  X10,9cm-2  for  curve  b.  I  Reprinted 
with  permission  from  L.M.Frontz  and  J.S. 
Nodvik,  J.  Appl.  Phys.34  ,  2348  (1963).] 


amplifier  dimensions.  Thus,  we  can  write  for  the  instantaneous  power16: 
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Fron.  the  power  expression,  we  can  estimate  the  expected  length  necessary  for  a  given  system 
and  an  effective  amplifier  cross  section.  It  can  be  shown  that  the  gain  and  energy  satisfy  the 
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TABLE  1-2 

SHORT-PULSE  POWER  AMPLIFIER  LENGTH 

A  (cm^) 

L  (m) 

0.01 

322 

0.03 

107 

0.10 

32 

0.30 

11 

1.00 

3.3 

3.00 

1.2 

10.00 

0.45 

30.00 

0.24 

100.00 

0.17 

We  have  used  a  program  implemented  by  S.  Marcus 
(see  Sec.  C-l-b)  to  calculate  the  inversion  for  a 
3:2:1  mixture  pumped  by  a  preionizing  e-beam  over 
a  range  of  primary  currents,  and  a  range  of  total 
pressure  holding  E/P  constant.  Table  1-2  gives 
estimates  of  the  amplifier  dimensions  one  might 
use  for  a  multistage  amplifier  system  to  bring  the 
individual  pulses  from  2  x  lo-"  to  2  J  total  energy 
("■nanosecond  pulses  and  10mA/cm2  primary  cur¬ 
rent;  the  length  restriction  i  somewhat  relaxed 
for  pulses  long  enough  to  tap  the  entire  energy 
reservoir  of  the  CO_-  vibrational  states).  In 
practice,  several  stages  would  be  used,  the  cross 
section  increasing  as  the  energy  in  the  pulse 
increased. 

At  this  juncture,  it  might  be  suggested  that  a 
solution  to  part  of  the  restrictions  is  to  run  the 
amplifier  system  unsaturated,  thus  eliminating  temporal  pulse  distortion  (the  unsaturated  re¬ 
quirement  being  more  easily  achieved  for  the  longer  FM  pulses).  If  the  duty  cycle  of  the  ampli¬ 
fier  is  small,  then  it  is  necessary  to  avoid  pun  ping  the  medium  at  hard  with  excitation  energy 
because  some  fraction  (~7S  percent)  will  always  go  into  medium  heating  and  the  remainder  will 
be  dissipated  by  this  means  also,  if  not  extracted  as  optical  energy.  Since  a  gas  (translational) 
temperature  rise  of  about  300"K  completely  shuts  off  the  gain  mechanism,  and  this  is  one  of  the 
two  significant  sources  of  pulse  distortion,  techniques  to  minimize  the  heating  over  the  range  of 
operational  modes  must  be  utilized.  For  the  long  coherence  time,  it  should  be  possible  to  re¬ 
duce  the  primary  pumping  (or  secondary  accelerating  field)  so  that  the  gain  buildup-to-optimum 
time  is  extended  to  about  the  time  that  the  next  pulse  arrives.  Thus,  medium  heating  is  reduced 
to  preserve  medium  quality  for  longer  periods.  Sufficiently  fast  flow  rates  would  aid  in  this  ef¬ 
fort.  Calculations  of  the  pumping  and  flow  requirements  are  presently  being  made.  A  serious 
source  of  medium  distortion  is  the  pressure  shock  wave  which  emanates  from  the  cathode.  One 
hypothesis  is  that  this  develops  from  sharp  temperature  differences  due  to  ion  acceleration  and 
then  ion-neutral  and  ion-cathode  collisions  in  the  cathode  fall  region.  Within  the  time  for  the 
wave  to  travel  at  the  speed  of  sound  transversely  across  the  amplifying  medium,  significant 
index-of- refract  ion  changes  can  occur  and  distort  the  medium.  Basic  research  is  needed  to 

study  this  effect  and  means  of  diminishing  or  eliminating  the  shock  during  the  operation  time  of 
the  amplifier. 

’a  tho  Flu  sequence,  another  potential  problem  that  was  investigated  was  the  artificial 

frequency  chirping  due  to  a  strongly  varying  index  of  refraction  (anomalous  dispersion)  in  the 
frequency  range  over  which  the  pulse  changes.  The  medium  index  of  the  CC>2  is  given  by44,4S 
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where  the  first  term  is  the  VHF  limit  of  t),  the  seeond  term  represents  the  significant  disper¬ 
sion  in  the  pain  region,  and  the  third  term  represents  the  contribution  from  all  other  higher  fre¬ 
quency  molecular  and  eleetronie  absorption  bands.  In  this  expression,  gR  is  the  line  center 
power  gain,  v  is  the  frequency  of  the  radiation,  r  is  the  line  center  gain  frequency,  Akq  is  the 
gain  half-bandwidth  (at  half-power),  and  e  is  the  speed  of  light.  The  seeond  term  could  con¬ 
ceivably  contribute  significantly  to  index  and,  hence,  wavenumber  changes  across  our  chirped 
pulse,  producing  a  phase  velocity  dependent  on  frequency,  and  thus  artificial  chirping.  Imagine 
that  we  send  in  a  pulse  whose  envelope  is  a  square  wave  of  duration  T  and  whose  fine  structure 
is  given  by 

I(t )  «  cos  { 2ir  [i-o  +  —■  (t  -  y)]  t) 


then  the  approximate  fine  structure  frequency  distribution  is  given  in  Fig.  1-48  vs  the  anomalous 
dispersion  of  the  medium  in  that  region. 


Fig.  1-48.  Frequency  content  of  chirped  pulse  in  region  of  resonant  gain  and  anomalous  dispersion. 

If  we  ask  the  change  in  ehirp,  we  find  that  a  chirp  of  6r/T  without  dispersion  becomes 

6 ii  fw  ,.  _  (it , 

T  +  .it  '  T  '  T' 

with  dispersion.  So  the  condition  for  the  change  in  ehirp  to  be  negligible  is  for 
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For  a  chirp  fir  =  2Ar  =  500  MHz,  and  a  length  of  1  m, 
ficfit  =  10"15 

which  is  indeed  negligible. 

Power  saturation  of  the  medium  is  seen  to  simply  smooth  out  the  local  dispersion  effect. 
Saturation  will  also  cause  some  bending  of  the  wavefront  at  a  fixed  transverse  cross  section  of 
the  pulse  if  the  intensity  strongly  varies  across  that  cross  section.  It  is  clear  that  a  power  dif¬ 
ference  twice  saturation  across  the  transverse  direction  will  produce  a  phase  shift  of  half  the 
magnitude  at  the  leading  edge  (fir  =  Ar)  across  the  beam.  The  index  change  due  to  dispersion, 
if  we  assume  fir  =  Ar  and  gQ  =  0.03  cm  *  =  13dB/m,  is  then 

~*go 

A^  +  Ar)  Sir- 
which  corresponds  to  a  phase  shift  of 

-Mg0z) 

2?rAnz  =  - ^ -  =  one  wavelength/17.  3  dB  of  gain 

(i.e.,  if  we  use  the  full  bandwidth  of  the  amplifier  for  chirping).  If  we  reduce  the  amplifier  band¬ 
width  used,  the  phase  shift  goes  down  less  than  linearly. 

H.  Granek 
R.  J.  Carbone 

2.  The  Possibility  of  Generating  Laser  Plasmas  by  Photoionization 

a.  Theoretical  Investigation 

The  possibility  of  producing  large  volume  uniform  plasmas,  via  photoionization,  with  elec- 

1 2 

tron  densities  on  the  order  of  10  /cc  and  greater  in  high-pressure  gas  laser  media  is  investi¬ 
gated  as  an  alternative  to  the  e-beam  ionization  technique.  As  in  the  e-beam-excited  lasers,  the 
electrons  generated  by  photoionization  would  be  heated  in  an  externally  applied  electric  field  to 
the  optimum  temperature  for  excitation  of  laser  states. 

Two  major  problems  are  encountered.  The  first  relates  to  the  magnitude  of  the  photoioniza- 

“48  2 

tion  cross  section  a •  For  most  atoms  and  molecules,  a  is  on  the  order  of  1  Mbn  =  10  cm 

(1  megabarn)  for  photon  energies  above  the  ionization  limit.  At  elevated  ionizing  gas  pressures, 

a  large  value  of  a  limits  the  incident  photons'  penetration  depth  to  a  short  distance.  For  instance, 

1  u  ,  3 

at  1  atm  the  gas  density  N  is  approximately  3X10  /cm  ,  and  for  (7^1  Mbn  the  penetration 
depth  L  =  i/Nct  is  «0.  3  mm.  In  this  case,  an  enormous  electron  density  is  produced  in  a  thin 
layer  of  the  gas,  which  may  be  of  use  in  special  applications  but  not  for  the  production  of  large- 
volume  plasmas. 

The  second  problem  to  he  considered  is  that  the  ionization  potentials  of  molecular  and  atomic 
species,  which  are  known  to  be  suitable  for  high-pressure  gas  lasers,  arc  generally  fairly  large. 

In  these  cases,  the  corresponding  photon  wavelengths  at  energies  above  the  ionization  limit  lie 
deep  in  the  vacuum  UV,  where  transparent  windows  are  not  readily  available  and  where  absorp¬ 
tion  processes  other  than  photoionizi.tion  may  become  important  in  nonionizing  species  in  the 
laser  medium. 
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Two  methods  are  proposed  to  overcome  these  problems.  The  first  is  to  seed  the  high- 
pressure  laser  medium  with  a  low-partial-pressure  gas  having  a  low-ionization  potential  to  allow 
ionization  by  photons  at  wavelengths  larger  than  UV  cutoff  limits  of  available  windows,  and  also 
longer  than  the  strong  UV  absorption  bands  of  the  high-pressure  host  laser  medium.  The  low 
partial  pressure  allows  deep  penetration  of  the  ionizing  photon,  and  hence  uniform  production  of 
photoelectrons.  An  example  of  this  approach  is  a  1-atm  COj-Ile-N^  laser  seeded  with  (1  1-torr 
cesium.  The  threshold  for  photoionization  of  cesium  is  3250A,  well  above  the  wavelengths  for 
the  onset  of  absorption  in  CO,  (1700A),  N,  (14S0A),  and  lie  (500A).  Calculation  shows  that, 
upon  irradiations  by  a  flashlamp  emitting  1  MW/cm  between  3250  A  and  2000  A ,  electron  density 
in  excess  of  t  0*  Vcm?  may  be  produced.  Seeding  with  cesium  may  also  be  applicable  in  the 
high-pressure  CO  laser.  By  generating  the  fourth  harmonic  of  a  Nd:  VAC.  laser,  it  is  possible 
to  produce  1  J  of  UV  in  10  nsec  at  2062  A.  This  wavelength  is  near  the  peak  of  the  ph  jtoionization 
cross  section  of  rubidium.  Using  this  intense  source  to  ionize  rubidium  may  produce  electron 
densities  of  S  X  10‘2/cni5.  There  are  a  host  of  organic  molecules  which  have  vapor  pressure  on 
the  order  of  1  torr  at  room  temperature  and  also  have  ionization  potentials  low  enough  so  that 
1'V  of  wavelengths  larger  than  1700A  may  produce  photoionization.  Among  these  are  tripentz- 
lamine,  tripropvlamine,  anti  tributylamine.  The  list  of  organic  compounds  with  low  ionization 
limits  is  so  extensive  that  a  search  is  expected  to  reveal  several  that  combine  the  necessary 
characteristics  of  high  vapor  pressure,  large  cross  section,  transparency  at  the  laser  wave¬ 
length,  and  low  ionization  potential. 

The  second  method  takes  advantage  of  the  inherently  large  photoionization  cross  section  of 
most  gases  and  attempts  to  induce  photoionization  by  means  of  a  two-step  (or  multiple-step)  exci¬ 
tation  where  two  (or  sc. oral)  incident  photons  are  absorbed  by  the  high-pressure  gas  to  produce 
a  single  photoelectron.  In  a  two-step  process,  the  photon  energies  can  be  as  low  as  about  half 
the  ionization  limit,  and  hence  the  window  transmission  problems  become  less  severe.  Further¬ 
more,  in  a  multiple-step  process,  the  reduced  size  of  the  effective  photoionization  cross  section 
would  allow  deep  penetration  of  the  incident  photons  and  production  of  electrons  over  a  largo 
volume.  This  method  may  be.  applicable  to  the  photoionization  of  one  of  the  molecular  or  atomic 
species  of  the  gas  laser  medium,  >r  it  may  he  applied  to  multiple-step  photoionization  of  an 
additional  gas  mixed  in  with  the  last  gas  laser  medium.  One  attractive  situation  is  found  in  CO 
where  the  lowest  lying  excited  electronic  states  are  triplets  lying  approximately  half-way  to  the 
ionization  limit.  The  transition  to  the  ground  state  is  only  weakly  allowed  so  the  absorption  cross 
section  is  small,  permitting  the  use  of  CO  at  1-atm  pressure.  The  intermediate  excited  states 

are  nearlv  metastable  with  a  long  radiative  lifetime  so  that  large  intermediate  state  densities  may 

4  () 

build  up  and  give  efficient  photoionization  in  the  second  step.  Calculation  shows  that,  with 
1-MW/cm2  illumination  between  1500A  and  2000 A,  electron  densities  on  the  order  of  10  /cm 

may  lie  produced  in  1  atm  of  CO.  j  <3  Levine 

A.  Javan 


b.  Experimental  Investigation 

An  ionization  cell  lias  been  constructed  that  contains  a  small  linear  Xe  flashlamp  and  two 
flat  parallel  electrodes.  The  lamp  is  oriented  parallel  to  the  electroties  in  their  midplane  and 
outside  of  the  uniform  field  region  which  exists  when  a  voltage  is  applied,  so  that  the  volume 
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between  the  electrodes  is  illuminated.  All  metal  corners  and  edges  are  removed  to  minimize 
ionization  due  to  corona  discharge.  The  cell  may  be  filled  with  various  buffer  gases  and  low 
ionization  potential  vapors  to  be  irradiated  with  an  intense  pulse  of  UV;  the  electrons  formed  by 
photoionization  are  drawn  off  by  the  electric  field  between  the  electrodes,  and  the  current  is 
measured  as  the  voltage  drop  across  a  sensing  resistor  in  series  with  the  high-voltage  supply. 
Because  the  drift  velocity  of  electrons  in  the  buffer  gas  (nitrogen)  is  known  as  a  function  of  elec¬ 
tric  field  divided  by  pressure  E/P,  the  electron  density  may  be  determined  from  the  current 
measurements  from  the  relation  j  =  neVj-y  where 

j  =  current  density  =  current/area  of  electrode 

e  =  electronic  charge 

v^j  =  drift  velocity  of  electrons  in  buffer  gas. 

The  Xe  lamp  was  supplied  by  the  Xenon  Corporation.  It  is  powered  by  a  Xenon  Corporation 
Micropulse r  and  can  be  fired  at  voltages  between  13  and  20kV;  the  maximum  voltage  represents 
an  input  energy  of  20  J.  The  lamp  arc  length  is  8  inches  (the  electrodes  arc  also  8  inches  long). 

O 

A  suprasil  envelope  was  used  to  permit  passage  of  UV  down  to  1650A,  although  it  is  expected  that 

O 

relatively  little  energy  is  emitted  at  that  wavelength;  the  output  peaks  around  2800  A.  The  light 
pulse  has  a  full  width  at  half-height  of  1  psec. 

When  the  cell  is  filled  with  N2  to  800  torr  and  a  potential  of  typically  1  0  kV  is  applied  to  the 

electrodes,  a  current  pulse  in  the  sensing  resistor  is  seen  which  closely  follows  the  light  pulse 

as  simultaneously  monitored  with  a  photomultiplier.  The  peak  current  corresponds  to  an  elec- 

tron  density  of  2  X  10  /cm  .  This  current  is  interpreted  as  the  result  of  photoemission  from  the 

"9 

aluminum  electrodes;  the  time  delay  of  less  than  10  sec  between  photon  absorption  and  electron 
emission  is  too  small  to  be  resolved. 

W'hen  0. 2-torr  tripropylamine  is  added  to  the  N.,,  a  dramatic  two-order-of-magnitude  in¬ 
crease  in  the  current  flow  is  seen  upon  flashlamp  ignition.  The  current  is  found  to  vary  linearly 
with  the  electrode  voltage  between  the  limits  corresponding  to  E/P  =  0  to  E/P  =  10.  This  implies 
that  a  peak  electron  density  on  the  order  of 
1  01  °/cm3  is  produced  independent  of  the  elec¬ 
tric  field.  In  turn,  this  is  strong  evidence  that 
the  enhancement  in  electron  production  is  due 
to  photoionization  and  not  to  ionization  of  tri¬ 
propylamine  by  electron  collision.  If  the  lat¬ 
ter  were  the  case,  then  an  increase  in  electron 
density  would  be  expected  as  E/P  (and  hence 
the  average  electron  energy)  is  increased; 
this,  in  fact,  is  not  observed. 

The  time  development  of  the  current  pulse 
is  markedly  different  from  the  case  of  pure 
N^.  As  seen  in  Fig.  1-49,  the  electron  density 
starts  to  grow  only  after  the  illumination  has 
reached  its  peak  and  continues  to  grow  until  the  light  pulse  terminates.  Thereafter,  the  electron 
density  decavs  -  rapidly  at  first,  and  then  with  a  long  exponential  tail.  These  observations  are 
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FLASHLAMP 


ELECTRON  DENSITY 


Fig.  1-49.  Upper  troce  shows  light  pulse 
Pom  floshlomp  meosured  by  1 P28  photo¬ 
multiplier  tube.  Lower  troce  is  current  flow 
between  electrodes  due  to  photoionizof ion. 
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tentatively  explained  as  follows.  Ionization  takes  place  by  a  two-step  process.  This  ionization 

O 

potential  is  about  7,5  V  so  that  two  2600-A  photons  have  sufficient  energy  to  liberate  an  electron. 
The  first  photon  populates  intermediate  states  half-way  to  the  ionization  limit;  the  second  photon 
then  raises  the  excited  electron  into  the  continuum.  The  buildup  of  intermediate  state  popula¬ 
tion  may  account  for  the  1-psec  delay  in  the  growth  of  the  electron  density.  Initial  rapid  decay 
of  the  plasma  maybe  due  to  recombination  which  proceeds  at  a  rate  proportional  to  n  .  When 
n  falls  to  a  sufficiently  low  value,  attachment  and  diffusion  become  the  dominant  loss  mechanisms 
which  proceed  at  a  rate  proportional  to  n  and  give  rise  to  an  exponential  decay. 

The  peak  electron  density  has  been  measured  as  a  function  of  flashlamp  power.  Preliminary 
results  show  that  for  relatively  low  lamp  energy,  the  electron  density  increases  approximately 
as  the  square  of  the  lamp  energy;  for  high  lamp  energy,  electron  density  increases  more  slowly  - 
approximately  linearly  with  lamp  energy.  This  behavior  is  consistent  with  two-step  excitation 
for  which  ionization  is  proportional  to  lamp  energy  squared.  At  low  lamp  energy,  and  hence 
low  electron  density,  electron  loss  varies  as  n;  thus,  the  equilibrium  electron  density  varies  as 

lamp  energy  squared.  At  high  lamp  energy,  and  hence  high  electron  density,  electron  loss  varies 
2 

as  n  ;  thus,  the  equilibrium  electron  density  varies  linearly  with  lamp  energy. 

Experimental  studies  are  being  pursued  to: 


(1)  Determine  quantitatively  the  electron  destruction  mechanisms  in 
tripropylamine, 

(2)  Study  the  photoionization  of  other  organic  molecules, 

(3)  Increase  plasma  densities  by  using  more  energetic  lamps, 

(4)  Produce  plasmas  via  photoionization  of  cesium  vapor, 

(5)  Produce  a  low-pressure  and  ultimately  an  atmospheric-pressure 
C02~N2-He  photoionization  excited  laser. 


J.  S.  Levine 


3.  Other  Device  Work 


a.  The  Feasibility  of  a  10-atm  CO,-N„-He  Laser 

L.  L, 

The  variation  in  10.6-pm  small-signal  gain  of  C02-N2-He  e-beam  excited  mixtures  as  a 
function  of  pressure  is  considered.  The  well-known  formula  for  small  signal  gain  is: 

-2  /  g. 


cv(v  =  V  ) 
o 


^^‘sponl^l/Z 


(N2  _  Nl) 


where 


a(v  =  vQ) 
N, 


N. 


spont 


Av 


1/2 


gain/centimeter  at  line  center 

population  density  in  upper  level:  g2  =  degeneracy 
population  density  in  lower  level:  g^  =  degeneracy 
spontaneous  transition  time  from  2  -*  1 
half-width  at  half-height  of  Lorentzian  gain  profile. 
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In  the  case  of  CW  pumping,  the  steady-state  population  density  hh  in  the  absence  of  a  satu¬ 
rating  optical  frequency  field  is  the  product  of  the  excitation  rate  lb  and  the  level  lifetime  r.. 
For  direct  electron  excitation  of  the  CO^  levels,  R.  is  proportional  to  the  electron  density,  the 
CO^  ground-state  density,  and  the  electron  collision  cross  section 


CO, 


R. 

1 


cm  v  N„ 
e  e  Gs 


A  similar  relation  holds  for  the  electronic  excitation  of  nitrogen  vibrational  levels.  In  both 
cases,  the  rate  of  excitation  is  directly  proportional  to  pressure.  The  rate  of  collisional  trans¬ 
fer  of  excitation  from  nitrogen  (r  =  1)  to  C02  (001)  is  proportional  to  the  excited-state  nitrogen 
density  N2  (p  =  1)  and  to  the  ground- state  CO^  density 

r  CO  CO 

R  ~N2(r  =  l)NGS2~N2(r  =  0)NGS2  . 


Since  N2  (v  =  1)  is  proportional  to  N2  (u  =  0),  the  collisional  transfer  rate  should  increase  with 
the  second  power  of  pressure.  At  lOatm,  the  10-pm  transitions  are  totally  pressure-broadened; 
so,  Ai^/2  ~  pressure.  Further,  level  lifetimes  are  inversely  proportional  to  pressure. 

So  far  it  appears  that,  in  the  CW  limit,  a  should  vary  inversely  with  pressure  with  respect 
to  direct  electron  excitation,  and  be  constant  in  pressure  with  respect  to  collisional  transfer  be¬ 
tween  N2  and  C02.  However,  in  the  case  of  e-beam  excitation,  ng  is  proportional  to  the  square 
root  of  the  energy  loss  of  the  ionizing  beam  to  the  gas  per  centimeter  (assuming  recombination 
is  the  dominant  mechanism  for  ion  loss).  This  stopping  power  is,  in  turn,  proportional  to  the 
gas  density  or  ng  ~  n/ pressure.  Therefore,  if  sufficient  primary  electron  energy  is  available 
to  penetrate  the  medium,  the  gain  of  a  high-pressure  (collisionally  broadened)  C02  laser  should 

exhibit  a  pressure  variation  that  is  the  sum  of  two  terms.  One  term,  due  to  direct  electron  ex- 

-1  /E 

citation,  varies  as  (pressure)  ;  the  other,  due  to  collisional  transfer  from  N,  to  CO,,  should 

+l/2  C  C. 

increase  as  (pressure)  . 

The  parameters  appearing  in  the  gain  equation  are  evaluated  for  an  8:1:1  mixture  of  He:CO,:N, 

4  2  E 

at  10-atm  pressure  assuming  longitudinal  e-beam  excitation  by  a  10  A,  3-nsec  e-beam  source 
(Febatron  706).  Peak  gain  for  the  P2Q  line  is  estimated  asE.E  percent/cm. 

J.  S.  Levine 

b.  Construction  of  a  10-atm  C02-N2-He  Laser 

A  laser  has  been  constructed  based  on  the  calculations  outlined  in  the  preceding  section.  Di¬ 
mensions  of  the  active  volume  are  1  x  1  x  20  cm.  Electron- beam  excitation  is  in  the  axial  direc¬ 
tion,  parallel  to  the  optical  axis.  The  e-beam  is  focused  by  means  of  a  magnetic  lens  placed 
between  the  Febatron  706  source  and  the  2-mil  titanium  electron  window  on  the  laser.  The  mir¬ 
ror  nearest  the  e-beam  source  is  formed  by  a  thin  aluminized  mylar  film  stretched  taut  over  an 
optically  flat  pyrex  annular  ring.  This  mirror  permits  free  passage  of  the  ionizing  beam  while 
providing  96-percent  reflectivity  at  10.6  pm.  The  far  end  of  the  laser  is  terminated  by  a  KC1 
Brewster  angle  window.  The  optical  cavity  is  completed  by  a  curved,  partially  transmitting  ger¬ 
manium  mirror. 

Further  development  of  this  system  awaits  the  availability  of  the  Febatron  source. 

•J.  S.  Levine 
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c.  Lamb  Dip  in  Sealed-Off  CO  leasers 

We  have  completed  measuring  and  evaluating  data  on  the  Lamb  dip  for  16  different  rotational- 
vibrational  transitions  in  sealed-off  CO  lasers  operating  at  a  total  pressure  of  1.5  torr.  The 
experimental  data  of  power  output  vs  frequency  were  matched  to  theoretical  curves  which  included 
relaxation  across  the  velocity  profile  as  a  major  factor  limiting  the  depth  of  the  Iamb  dip.  As  a 
result  of  these  investigations,  the  following  parameters  were  accurately  determined:  the  inhomo¬ 
geneous  and  homogeneous  linewidths,  the  relaxation  rate  out  of  the  lasing  level,  and  the  cross¬ 
relaxation  rate.  We  found  that  the  cross-relaxation  rates  were  of  the  same  order  in  magnitude 
as  the  dephasing  time.  A  full  report  detailing  our  work  is  nearing  completion. 

C.  Freer! 

II.  A  Haus 
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II.  OPTICAL  MEASUREMENTS  AND  INSTRUMENTATION 


A.  ZnS  UP-CONVERSION 

Additional  parameters  relevant  to  the  up-conversion  process  in  ZnS  have  been  measured. 
Whereas  previous  measurements  have  been  primarily  concerned  with  the  effects  of  DC  IR  radia¬ 
tion,  AC  effects  have  also  heen  measured  now  using  chopped  IR  radiation  as  a  source  and  con¬ 
stant  UV  illumination  for  excitation.  Under  these  conditions,  an  AC  visible  output  is  observed 
on  a  DC  background.  At  frequencies  on  the  order  of  typical  TV  scanning  rates,  the  following 
results  were  obtained: 

(t )  Peak-to-peak  AC  visible  output  decreases  about  50  percent  in  going 
from  1  to  64  Hz; 

(2)  AC  visible  output  increases  with  increasing  UV  intensity; 

(3)  AC  to  DC  background  ratio  decreases  with  increasing  UV  intensity; 

(4)  AC  signal-to-noise  ratio  (SNR)  increases  somewhat  with  increasing 
UV  intensity; 

(5)  AC  SNR  is  1  to  2  orders  of  magnitude  worse  than  the  SNR  obtained  in 
previous  measurements  where  the  phosphor  is  exposed  to  DC  IR  radia¬ 
tion  after  (not  during)  UV  excitation.  This  is  due  to  the  large  back¬ 
ground  noise  caused  by  constant  UV  illumination  in  the  AC  case. 

These  results,  and  others,  simply  confirm  the  fact  that  in  these  materials  it  takes  a  finite  time 
(several  seconds)  to  excite  electrons  to  states  which  will  produce  the  maximum  stimulation 
effect.  Under  AC  operation,  even  at  low  frequencies,  there  is  not  enough  time  for  significant 
UV  excitation  between  IR  pulses.  In  other  measurements,  we  have  determined  that  the  spectral 
character  of  the  visible  output  is  dependent  upon  the  IR  stimulation  wavelength.  Our  results 
show  that  10-pm  radiation  produces  a  band  of  visible  radiation  centered  near  5300  A,  while  2-pm 
radiation  produces  a  band  centered  nearer  5400  A. 

Presently,  we  are  constructing  apparatus  to  determine  the  S/N  characteristics  of  the  phos¬ 
phor  screen  resolution  elements  under  uniform  IR  illumination.  We  expect  to  determine  the 
degree  of  nonuniformity  in  response  due  to  irregularities  in  particle  size,  particle  distribution 
and  layer  thickness.  The  smallest  particle  size  range  is  1  0  to  20pm. 

W.  J.  Scouler 

B.  IR  CAMERA  (4  inch) 

In  measurements  relevant  to  the  imaging  quality  of  the  ZnS  up-converter  phosphors,  we 
have  used  the  cooled,  4-inch-diameter,  f/l.O,  3-element-refractor  camera  shown  schematically 
in  Fig.  II - 1 . 

The  camera  has  a  0.25-rad  FOV  and  is  diffraction  limited  over  the  design  wavelength  range 
of  7.5  to  12  pm.*  A  refractor  was  chosen  over  a  reflector  for  the  f/l.O  design  as  it  has  a  simpler 
and  more  compact  configuration. 


*  Optic.il  design  In  l).i\id  (Jrr\  Associates,  Waltham,  Mass. 


Preceding  page  blank 
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PHOSPHOR  ON 
SAPPHIRE  SUBSTRATE 
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Fig.ll-1.  Schematic  of  4-inch  f/1 ,  0.25-rad  FOV,  diffraction  limited  I R  camera. 


60 


Section  II 


An  important  consideration  in  the  design  is  the  reduction  to  a  minimum  of  extraneous  lit 
radiation  at  the  image  plane.  This  was  achieved  by  cooling  I,j  to  liquid-lie  temperature  and 
to  liquid-N^  temperature-  and  then  mounting  between  them  a  shutter  cooled  to  nearly  1  iquid-IIe 
temperature  and  activated  via  an  O-ring  seal  at  room  temperature.  Cooling  to  liquid-N^ 
temperature  reduces  the  heat  load  to  the  liquid-lie  section  while  reducing  emission  from  I 
The  cooled  shutter  prevents  any  front-end  IR  radiation  from  reaching  the  phosphor  until  the 
scene  is  to  be  viewed. 

Behind  the  phosphor  screen  are  IR  blocking  filters  which  absorb  radiation  greater  than  2  pm 
while  being  highly  transmitting  to  both  visible  and  UV  radiation  Each  set  of  filters,  one  at 
liquid-N.,  and  one  at  liquid-lie  temperature,  consists  of  one  3-mm-thick  Oriel  7100  filter  and 

<-  -5 

one  3-mm-thick  piece  of  plexiglas.  The  attenuation  of  IR  2  pm  is  about  10  in  each  set.  The 
UV  excitation  light  and  the  visible  luminescence  from  the  phosphor  are  viewed  via  the  reflex 
camera  used  to  photograph  the  phosphor  image. 

The  4-inch-diameter  Ge  collector  lens  L.  also  serves  as  the  vacuum  window.  A  simple 
stress  analysis  shows  that  the  atmospheric  pressure  produces  no  measurahle  deterioration  of 
the  diffraction  limited  image.  L.  is  held  in  an  adjustable  holder  to  permit  some  focusing  ability 
for  objects  closer  than  infinity.  A  bandpass  filter  directly  in  front  of  the  phosphor  can  be  used 
to  limit  the  wavelength  range  to  any  band  hetween  7.5  and  12  pm. 

The  reservoirs  for  the  coolants  are  1.9  and  2.7  liters,  respectively,  for  liquid  He  and  liquid 
N_  giving  an  operating  time  in  excess  of  4j  hours.  Typical  cool  down  times  are  on  the  order  of 
3  hours. 

Mechanical  rigidity  and  thermal  isolation  of  the  three  different  temperature  sections  are 
achieved  using  specially  fabricated  ribhed  stainless  steel  tubing.  These  isolators  have  a  high 
resistance  to  twist  or  distortion  and  a  high  thermal  resistance. 

Tests  have  demonstrated  that  the  mechanical  integrity  of  the  design  is  maintained  at  cryo¬ 
genic  temperatures  and  that  its  performance  is  diffraction  limited.  This  was  determined  exper¬ 
imentally  hy  replacing  the  ZnS  phosphor  screen  with  a  Ge:Cu  detector  located  behind  a  knife 
edge  in  the  focal  plane,  looking  at  a  chopped  hlackbody  source  at  different  pointing  angles.  The 
variation  in  angle  from  cut-on  to  maximum  signal  is  a  measure  of  the  image  size.  From  our 
measurements  of  a  1-rm  source  at  45  m,  the  best  image  obtained  had  a  38-pm  diameter.  This 
is  to  be  compared  to  a  theoretically  calculated  value  of  34  pm,  which  is  within  experimental  limits. 

Rotation  of  the  camera  through  large  angles  produces  the  signal  vs  angle  curve  of  Fig,  II- 2. 
This  curve  was  obtained  with  a  Ge:('u  detector  located  behind  a  mask  at  the  focal  plane  defining 
an  area  on  the  detector  of  1.5  x  3  mm.  The  width  of  the  central  peak  corresponds  to  the  1.5-mm 
slit  width,  which  is  much  larger  than  the  image  size.  Thus,  the  detector  acceptance  area  for 
scattered  radiation  is  much  larger  than  is  necessary  to  view  the  direct  image.  If  the  detector 
area  were  reduced  to  just  accept  the  direct  image,  then  the  off-axis  signal  would  he  about  3  orders 
of  magnitude  lower  than  shown.  Note  that  at  angles  of  ahout  18"  from  the  normal,  there  are 
peaks  at  about  1/3000  of  the  signal  detected  with  the  camera  pointed  at  the  source.  Some  effort 
has  been  made  experimentally  (so  far  unsuccessfully)  to  determine  the  cause  of  these  peaks.  In 
addition,  computer  calculated  ray-traces  of  several  assumed  reflection  possibilities  have  pro¬ 
duced  no  equivalent  path  of  suspec  ;od  rays. 
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Fig.  11-2.  Detector  signol  vs  ongle.  Angle  in  degrees  is  omount  blockbody  source  is  located 
off  optical  axis. 


Presently,  a  second  set  of  lenses  ( ,  L 2  and  L^)  is  being  set  up  on  an  optical  bench  at 
room  temperature  in  an  attempt  to  discover  the  path  of  the  extraneous  rays. 

D.  If.  Dickey 
T.  M.  Quist 
W.  J.  Scouler 


C.  INTERFEROMETRIC  IMAGING  TECHNIQUES 
1.  Introduction 

The  most  serious  limitation  to  using  the  interferometric  imaging  technique  from  the  KC-135 
aircraft  in  flight  is  suspected  to  be  the  presence  of  an  image  degrading  boundary  layer.  The 
earlier  reported  retroreflector  type  measurements  showed  that  an  image  degradation  arises  be¬ 
tween  the  wing  tip  of  the  aircraft  and  the  skin.^ 

In  these  experiments,  the  measured  modulation  transfer  function  (MTF)  fell  quite  rapidly  to 
low  contrast  values,  and  there  was  some  question  raised  as  to  whether  the  MTF  actually  repre¬ 
sented  only  the  boundary  layer  (near  the  skin  of  the  aircraft).  In  order  to  confirm  that  the  degra¬ 
dation  was  in  fact  localized  near  to  and  just  outside  the  aircraft  window,  a  new  experiment  was 
designed  and  carried  out.* 

The  boundary  layer  localized  in  the  first  25  cm  away  from  the  skin  of  the  aircraft  was  meas¬ 
ured  in  this  experiment.  In  so  doing,  it  was  shown  that  the  interferometer  performed  properly 

*  I'Apt'riliKMtl  |>[.  niiicil  .mil  implemented  in  coopcrution  with  I).  korlieroml  II.  tt.  Bn  mil  in  support  of  the  I.KI.d  telescope  seeinp 
ussessineiil  work,  lint  usiiip  the  interferometer  os  o  sensor. 
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out  to  the  highest  spatial  frequency  (maximum  shear  value)  and  demonstrated  that  the  effects  of 
vibration,  air  currents  inside  the  aircraft  cabin  near  the  window  and  temperature  gradients 
across  the  optical  window  produced  no  significant  part  of  this  degradation. 

MTF  measurements  of  the  boundary  layer  were  made  both  at  the  front  (LFLC)  window  and 
at  the  third  window  in  the  aft  array.  The  measurements  were  then  compared  with  the  earlier 
retroreflector  results. 

2.  Experimental  Program 

A  schematic  of  the  experimental  setup  is  shown  in  Fig.  II- 3.  A  5-mW  He-Ne  CW  laser 
(wavelength  632. 8nm)  was  mounted  rigidly  on  top  of  the  interferometer.  A  high  quality  50-mm- 
aperture  laser  beam  expanding  telescope  was  used  witli  the  laser  to  provide  a  collimated  and 
coherent  laser  beam.  A  mirror  and  a  beam  splitter,  mounted  to  allow  angular  adjustments 
in  each  case,  were  arranged  to  direct  the  50-mm  collimated  laser  beam  through  a  75-mm- 
diameter  optical  window  in  the  side  of  the  aircraft.  The  beam  then  propagated  through  the 


Fig.  11-3.  Schematic  af  interferometric  boundary  loyer  MTF  measurement. 

nirstrenm  (the  boundary  layer)  to  a  gold-coated  mirror  31  mm  in  diameter  (supported  on  an 
airfoil  outside  the  aircraft)  where  it  was  reflected  back  through  the  boundary  layer  and  the  optical 
window.  The  mirro-  ’-as  located  25  cm  from  the  window  face  outside  the  aircraft.  Both  the 
window  and  mirror  wi  re  mounted  in  an  aerodynamically  designed  airfoil  plate  unit  made  of  alumi¬ 
num,  which  can  be  fitted  into  any  of  the  aircraft  window  locations  in  place  of  the  standard  22- inch- 
diameter  optical  windows  of  the  KC-135  aircraft.*  An  aerodynamically  designed  tube  (an  airflow 

*  Xil'f'iil  |d.it<*  design  In  11.  Mile*.  \irlit<nit‘  KcsimitIi  \  !\n^inrrmi^.  I!.in*ium  M'  ll. 
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shield)  could  also  be  readily  fitted  (on  the  ground)  between  the  gold-coated  mirror  and  the  75-mtr 
optical  window  to  completely  enclose  and  shield  the  light  path  outside  the  aircraft  from  the  bound¬ 
ary  layer. 

The  laser  beam  returned  through  the  optical  window,  following  the  double  pass  through  the 
boundary  layer,  and  passed  next  through  the  beam  splitter  to  a  second  smaller  24-mm-aperture 
beam  reducing  telescope.  The  beam  was  reduced  in  diameter  to  about  an  8-mm  collimated  beam 
and  then  passed  through  a  dove  prism  and  finally  to  the  corner  cube  shearing  interferometer. 

The  dove  prism  rotated  the  wavefront  around  the  optical  axis,  to  allow  measurements  for  differ¬ 
ent  diametrical  orientations  of  the  shear  axis.  The  diameter  of  the  collimated  laser  beam 
through  the  boundary  layer  was  limited  by  the  diameter  of  the  entrance  aperture  of  the  telescope 


to  a  nominal  24  mm. 

All  the  auxiliary  optical  components  arranged  in  this  experimental  setup  were  of  very  high 
optical  quality,  to  ensure  that  the  wavefront  entering  the  interferometer,  in  the  absence  of  the 


boundary  layer,  was  plane  to  within  l/lO  wave. 

A  photograph  of  the  interferometer  and  the  laser  source,  with  the  auxiliary  optical  compo¬ 
nents  in  place,  is  shown  in  Fig.  ll-4(a).  The  complete  unit  was  mounted  on  a  rigid  tripod,  which 
allowed  angular  and  translational  adjustments  to  be  made  and  allowed  the  laser  beam  to  be  di¬ 
rected  through  the  75-mm-diameter  window  in  the  airfoil  plate  (a  paper  mask  in  the  photograph 
shows  the  location  of  the  airfoil  plate  and  the  two  windows).  A  photograph  of  the  airfoil  plate  in 
position  outside  the  aircraft  is  shown  in  Fig.  ll-4(b).  The  laser  beam  passes  and  returns  through 
the  left  hand  window  in  Fig.  ll-4(b).  The  tripod  was  bolted  to  the  floor  *o  produce  a  relatively 
rigid  mounting.  The  laser  beam  was  aimed  through  the  right  hand  window  as  seen  in  Fig.  ll-4(a) 
and  a  second  window  (to  the  left  of  the  laser  beam  window)  in  the  airfoil  plate  served  simply  for 


visual  observation. 


Fig.  II -4(a).  Interferometer  at  LFLC  window 
position  in  aircraft. 


Fig.  I  l-4(b).  Airfoil  plate  from  outside 
oircroft  (airflow  shield  removed). 
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All  the  electrical  measurements  were  recorded  on  an  FM  tape  recorder  and  the  data  were 
analyzed  on  the  ground  in  the  laboratory.  An  oscilloscope  display  of  all  the  signals  provided 
continuous  real-time  monitoring  of  the  data  being  recorded  and  permitted  accurate  electro-optic 
focusing  of  the  telescope  in  flight  before  each  MTF  run. 

Both  of  the  two  airfoil  plate  windows  had  provision  for  either  a  single  or  double  (airspaced) 
optical  window  mounting  (to  circumvent  frosting  problems),  but  only  the  single  window  unit  was 
employed  in  the  experiment  reported  here. 

3.  Results 

Seven  daylight  flights  (at  about  40,000  feet  altitude)  were  carried  out.  In-flight  MTF  meas¬ 
urements  were  performed  with  the  airflow  shield  in  place  (no  boundary  layer)  and  with  it  removed 
(boundary-layer  measurement).  An  earlier  check  of  the  complete  system  was  done  in  the  labora¬ 
tory  and  in  the  aircraft  on  the  ground.  MTF  measurements  in  flight  were  made  at  the  LFLC 
window  position  (located  well  forward  of  the  wing,  7  m  from  the  aircraft  nose)  and  at  the  third 
mount  aft  array  window  position  (located  over  the  wing,  18  m  from  the  aircraft  nose).  The  inter¬ 
ferometer  was  tripod-mounted  for  the  LFLC  window  measurement,  and  was  placed  on  the  fixed 
gimballed  mount  (with  micrometer  adjustments)  for  the  third  mount  measurement. 

MTF  measurements  were  made  for  different  orientations  (2<p)  of  the  wavefront,  correspond¬ 
ing  to  the  MTF  along  a  horizontal  (2 <p  ~  0),  inclined  (2 <p  =  45”  )  and  vertical  (2 <p  =  90° )  axis.  This 
will  suffice  to  establish  whether  there  is  any  asymmetry  in  the  wavefront  arising  from  a  non- 
symmetry  in  the  boundary  layer.  Bach  telescope  could  be  independently  and  accurately  focused, 
and  the  focus  was  checked  before  each  MTF  run  by  maximizing  the  signal  on  the  scope  display 
at  a  shear  setting  of  approximately  s  l.n.  The  parameter  s  is  defined  as  a  reduced  spatial 
frequency  s  =  2XRF,  where  X  is  the  wavelength  (in  cm),  R  is  a  spatial  frequency  (in  cycles 
per  cm),  and  F  is  the  ratio  of  the  focal  length  to  aperture  diameter  of  the  optics. 

The  MTF  measurements  obtained  at  the  third  mount  aft  array  window  are  shown  in  Fig.  11-5 
(from  3  flights).  A  similar  set  of  data  were  obtained  at  the  LFLC  position,  and  the  respective 
mean  MTF  for  each  window  position,  is  shown  compared  in  Fig.  11-6.  In  both  figures  the  diffrac¬ 
tion  limited  (clear  aperture)  MTF  curve  is  shown.  The  data  represent  a  double  pass  through 
the  boundary  layer. 

With  the  airflow  tube  in  place  (no  boundary  layer),  the  MTF  is  close  to  the  diffraction  lim¬ 
ited  curve  (curve  1,  Fig.  11-5).  When  the  tube  is  removed  (boundary  layer  present)  the  MTF  falls 
(curve  2,  Fig.  11-5).  These  data  characterize  the  image  degradation  due  to  the  boundary  layer 
during  daylight  flights  at  40,000  feet  altitude  and  at  normal  cruising  speed.  There  is  no  signifi¬ 
cant  difference  in  data  obtained  for  different  pupil  rotations  (at  0”,  45°  and  90”)  corresponding  to 
different  shear  orientations  along  the  wavefront  diameter.  There  appears  therefore  to  be  no 
asymmetry  in  the  degradation  due  to  the  boundary  layer. 

Since  with  the  air  flow  tube  in  place  the  MTF  measured  is  virtually  diffraction  limited,  tnere 
cannot  be  any  significant  degradation  aF  ing  from  extraneous  air  currents  near  the  window  inside 
the  aircraft,  nor  from  temperature  gradients  in  the  optical  window  itself. 

In  order  to  compare  these  data  with  the  theoretical  prediction  available  at  this  time,  the 
characteristics  of  the  refractive  index  (or  temperature)  fluctuations  in  the  turbulent  airstream 
must  be  evaluated.  These  in  turn  will  determine  the  random  wavefront  errors  from  which  the 
MTF  can  be  computet!  for  this  case. 
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Fig.  11-5.  MTF  af  aircraft  boundary  layer  at  third  Fig.  11-6.  MTF  af  in-f I  ight  c  ircraft  boundary  layer 

windaw/aft  array.  at  front  (LFLC)  and  third  windaw/aft  array. 

As  a  rough  estimate  of  the  wavefront  errors  in  the  present  instance  for  a  double  pass  through 
the  aircraft  boundarv  layer,  a  fluctuation  for  the  optical  phase  (path  length)  of  X/7. 5  is  believed 
to  be  of  the  right  order*  This  value  will  be  used  without  detailed  discussion  to  make  a  prelim¬ 
inary  comparison  with  theory.  Barakat3  has  shown  that  the  MTF  can  be  calculated  from  random 
wavefront  errors.  Using  his  method,  the  correlation  parameter4  was  roughly  estimated  from 
our  earlier  retroreflector  MTF  measurements.  For  v  -  X/7.  5,  the  MTF  f^(s)]  was  then  computed 
from  the  relationship  given  by  Barakat 

VTsj  ~  exp{_k2ff2  t1  “  exp  (-as2)]} 

where  ^0<s)  is  the  diffraction-limited  MTF,  s  is  the  reduced  spatial  frequency  and  k  =  2tt/x. 

This  computed  curve,  for  the  conditions  spelled  out  above,  is  shown  by  the  points  on  Fig.  II-6. 
While  the  present  analysis  is  incomplete  and  should  be  regarded  as  tentative,  a  comparison  be¬ 
tween  the  measured  and  computed  curve  shows  satisfactory  agreement. 

In  order  to  compare  the  present  results  with  those  previously  obtained  in  the  wing-tip  retro- 
reflcctor  experiment,  in  which  data  for  effectively  a  single  pass  were  obtained  with  an  89-mm- 
diameter  heam  (with  a  central  obscuration  «  =  0.34),  a  certain  amount  of  extrapolation  is  called 
for. 

The  equivalent  single-pass  MTF  through  the  boundary  layer,  can  be  obtained  from  the  actual 
measured  double-pass  data  (Fig.  II- 5 )  by  assuming  that  the' relationship  between  the  MTF  (ip)  and 
the  path  length  L  through  the  boundary  layer  is  of  the  form:  L  =  k  log (^Q/$)  where  k  is  a  con¬ 
stant  and  ?/'  is  the  optical  diffraction  limited  MTF.  The  equivalent  single- pass  MTF,  computed 


*  l.stim.iled  frimi  sheared  interferngrams  photographed  during  the  airi  raft  flights. 
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in  this  way  from  the  double-pass  measured  curve,  is  next  extrapolated  from  the  24-mm-diameter 
clear  aperture  case  to  the  89-mm-diameter  obscured  aperture  case. 

In  this  extrapolation,  it  is  assumed  that  there  is  no  inherent  aperture  size  dependence,  except 
for  scaling,  on  the  degradation.  It  is  simply  an  extrapolation  along  the  spatial  frequency  axis,  to 
correlate  this  parameter  with  the  respective  beam  diameter.  The  resulting  curve  from  this 
extrapolation  is  shown  in  Fig.  II- 7,  in  the  case  of  the  third  window/aft  array  measurement.  It  is 
compared  on  the  same  figure  with  the  MTF  measurement  made  in  the  retroreflector  experiment 
at  this  same  window  (daytime  flying).  For  the  lower  spatial  frequencies  (below  s  =  0.2)  there  is 
relative  agreement  between  the  two  curves.  The  extrapolated  curve  above  s  =  0. 2  is  not  very 


Fig.  11-7.  MTF  of  in-flight  oircroft  boundory 
layer  ot  third  window/oft  orroy  [comparison  of 
direct  (25-cm)  path  ond  retroreflector  (19-m) 


accurate,  since  it  is  determined  from  data  at  the  higher  spatial  frequencies  in  the  24-mm-beam 
experiment  where  the  SNR  is  poor.  Nevertheless,  it  provides  a  direct  comparison  of  the  in¬ 
flight  MTF  measurement  for  the  local  (25-cm)  path  and  the  retroreflector  experiment  (19m  to 
the  wing  tip)  path  at  the  third  window/aft  array  location.  A  similar  comparison  was  made  for 
the  LFLC  position,  with  slightly  poorer  agreement  between  the  two  curves,  due  mainly  to  the  in¬ 
fluence  of  window  astigmatism  on  the  retroreflector  MTF  curve  at  the  LFLC  position. 

4.  Conclusions 

The  boundary  layer  is  mainly  responsible  for  the  degraded  imaging  observed  in  flight  through 
the  two  aircraft  windows  mentioned,  the  seeing  being  slightly  better  at  the  front  window.  By  com¬ 
paring  the  local  (25-cm)  path  measurements  with  the  longer  (19-  or  27-m  retroreflector)  path 
measurements  it  is  concluded  that  the  degradation  apparently  is  confined  mostly  within  the  25-cm 
path. 

An  attempt  made  to  obtain  theoretically  computed  MTF  curves,  based  upon  present  knowledge 
(this  ultimately  has  come  from  wind  tunnel  data  to  date),  gives  reasonable  agreement  between 
theory  and  experiment  for  the  airplane.  The  assumptions  and  estimations  made  here  should  be 
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regarded  as  tentative,  until  better  substantiation  of  the  value  for  the  wavefront  errors  is  obtained. 

In  summary,  the  interferometric  system  has  been  demonstrated  to  be  capable  of  MTF  meas¬ 
urements  over  its  complete  range  (to  maximum  shear)  on  board  an  aircraft  in  flight,  and  has 
established  the  existence  of  image  degradation  due  to  a  boundary  layer. 

D.  Kelsall 

D.  AIRBORNE  IMAGE  STUDIES 

1.  Measurements  of  the  Tropopause  Effect 

In  the  fall  of  1971,  four  flights  were  made  with  the  Lincoln  Laboratory  KC-135  in  the  neighbor¬ 
hood  of  Christchurch,  New  Zealand,  in  order  to  obtain  high  resolution  stellar  photographs  both 
above  and  below  the  atmospheric  tropopause.  The  purpose  was  to  determine  whether  any  effect 
could  be  observed  on  the  images  caused  by  the  turbulent  air  in  the  boundary  region  and  to  relate 
the  measurements,  if  possible,  to  proposed  profiles.  D.  Kelsall  has  reported^  the  results  of 
MTF  measurements  which  were  made  on  the  same  flights.  He  noted  no  significant  trend  in  the 
MTF  above  and  below  the  tropopause,  nor  any  elevation  angle  dependence  (Ref.  5,  p.  47). 

A  30-cm-diameter  Schmidt-Cassegrainian  diffraction-limited  telescope  was  used  to  photo¬ 
graph  the  images  on  Kodak  plus-X  film.  Images  formed  by  this  telescope  in  flight  have  been  ex¬ 
tensively  analyzed;  the  various  sources  of  image  degradation  have  been  identified  and  measured. 

In  this  series  of  experiments  we  had  hoped  to  extract  the  "tropopause  component"  T  in  the  same 
way  that  the  other  components  were  found,  namely  by  root  differences  of  squares.  Specifically, 
two  separate  experiments  were  performed.  The  first  consisted  simply  of  photographing  a  star 
from  5  kft  below  the  tropopause  and  then  from  5  kft  above;  a  complete  focus  run  was  made  at  each 
altitude.  Determining  the  average  altitude  of  the  tropopause  was  facilitated  by  balloon  soundings 
received  every  12  hours  from  several  locations  within  a  radius  of  100  to  200  miles.  The  Christ¬ 
church  meteorological  office  was  particularly  helpful  in  obtaining  the  necessary  observations  and 
forecasting  information.  The  estimated  tropopause  position  was  checked  by  recording  the  free 
stream  temperature  vs  altitude  on  the  aircraft  during  the  experiment  and  by  noting  the  altitude  at 
which  temperature  inversion  occurred.  The  second  experiment  consisted  of  photographing  a  star 
for  10  sec  every  minute  as  the  aircraft  climbed  through  t.he  "disturbed"  region.  This  was  re¬ 
peated  at  a  predetermined  focus  setting  on  each  flight.  The  rate  of  climb  varied,  but  averaged 
about  1000  feet/min.;  the  total  climb  was  about  10,000  feet.  The  main  advantage  of  the  second 
experiment  over  the  first  was  that  local  weather  changes  could  be  kept  to  a  minimum  during  pas¬ 
sage  through  the  tropopause. 

After  the  photographs  were  developed,  about  5  images  were  selected  at  random  under  each 
condition,  enlarged  approximately  400X,  and  their  diameters  measured  by  estimating  the  diam¬ 
eter  containing  9  0  percent  of  the  exposed  grains  using  calibrated  circles.  This  method  has  been 
calibrated  by  an  actual  grain  count  to  establish  the  "eye  ball"  9  0-percent-circle.  Two  sets  of 
independent  measurements  were  made,  one  by  D.  G.  Kocher  and  the  other  by  B.  W.  Bryant.  The 
results  were  averaged  to  minimize  observer  bias.  The  results  of  the  4  flights  are  shown  in 
Fig.  II- 8.  There  appears  to  be  little  if  any  dependence  of  image  size  on  the  tropopause  position 
for  any  of  the  flights.  Other  disturbances,  possibly  local  weather,  account  for  most  of  the 
changes  in  image  size.  For  example,  the  unusually  fine  images  on  flight  3  at  about  14°  elevation 
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Fig.  11-8.  Image  size  vs  elevation  angle,  above  and  below  the  tropopause.  Flights  out 
of  Christchurch,  New  Zealand,  September  1971. 


are  otherwise  very  difficult  to  understand.  Of  course,  the  tropopause  is  a  very  ill-defined 
region,  in  general;  the  strength  of  the  atmospheric  turbulence  there  can  vary  greatly.  Conse¬ 
quently,  one  would  certainly  not  expect  to  find  a  constant  tropopause  component.  If  we  let 
be  image  size  below  the  tropopause  and  <z>A  be  that  above,  then  the  mean  value  of  the  component 
T  -  -  (Ac?.2)1/2.  However,  from  Fig.  11-8,  one  finds  that  A <p}  is  negative  (T  imag¬ 

inary)  about  40  percent  of  the  time.  Only  at  the  largest  elevation  angles  do  there  seem  to  be 
consistently  better  images  above  the  tropopause  with  T  =*  7  X  10  ^  rad,  but  this  is  probably  il-  ^ 
lusory  since  7  x  10"&  is  about  the  lower  limit  of  statistical  significance.  In  fact,  if  there  were 
a  significant  T-component  on  any  of  our  flights,  one  would  have  expected  to  see  a  greater  and 
more  consistent  difference  at  low  elevation  angles  where  the  path  through  the  disturbed  region 
is  greater.  If  one  is  located  just  at  the  base  of  the  disturbed  region,  then  the  path  length  varies 
approximately  as  cosec  of  the  elevation  angle.  The  path  at  4°  elevation  is  about  1  OX  what  it  is 
at  50°  and  therefore  T(4° )  should  be  about  3X  T(50°  ).  All  of  our  climb-out  data  were  obtained 
at  rather  low  elevations  in  order  to  gain  as  much  sensitivity  to  the  turbulent  region  as  possible. 
Figures  11-9  through  -12  summarize  these  data.  Again,  the  fluctuations  appear  to  be  independent 
of  the  position  of  the  tropopause  within  the  sensitivity  of  the  measurements.  Only  in  flight  4 
(Fig.  11-12)  where  the  elevation  angle  was  4.1°  is  there  even  a  suggestion  of  smaller  images  in 
the  stratosphere.  These  data  lead  to  T  ^  13  X  10  ^-rad  as  an  upper  limit  for  the  contribution  of 
the  stratospheric  boundary- layer  turbulence  to  the  90-ptu  cent  image  size. 

Admittedly,  these  data  represent  a  small  sample  at  one  location  and  time  of  year.  However, 
similar  data  were  Collected  near  Anchorage,  Alaska,  in  August  1963.  Though  not  as  extensive, 
the  results  were  very  similar  to  those  shown  in  Fig,  11-8,  The  main  limitation  to  the  sensitivity 
of  this  method  is  the  existence  of  a  large  and  variable  "airborne  seeing"  component  (now  identified 
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Fig.  11-9.  Image  size  vs  altitude,  flight  1, 
4  September  1971 


Fio ,  11-11.  Image  size  vs  altitude,  flight  3, 
9  September  1 971 . 
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Fig.  11-10.  Image  size  vs  altitude,  flight  2, 
6  September  1971 . 
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12  September  1971 . 
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principally  with  the  aircraft  boundary  layer),  which  causes  the  images  to  be  much  larger 
(~20  x  10  rad)  than  the  Airy  disk  (~4  x  lo"6  rad).  The  causes  for  its  fluctuation  are  not  known, 
but  may  be  related  to  the  air  mass  changes  in  the  upper  atmosphere.  During  relatively  quiet 
periods,  the  sensitivity  limit  is  determined  by  the  statistical  fluctuations  in  the  measurements 

1  x  1 0  rad )  of  these  "enlarged"  images.  Therefore,  the  limiting  sensitivity,  Sj.mit  = 

10  v(21)  -(20)  =6.3X10  rad;  i.e.,  any  component  less  than  this  is  lost  in  the  "noise." 

Because  additional  noise  sources  are  always  present,  the  practical  limit  is  probably  closer  to 
8  or  9  x  10  rad.  These  results  are  currently  being  compared  with  various  theoretical  models 
of  the  upper  atr.oaplicre. 

2.  Boundary  Layer 

Another  series  of  aircraft  experiments  was  planned  and  carried  out  in  March  and  April  of 
1972  which  showed  convincingly  that  the  turbulent  boundary  layer,  a  region  about  10  to  25  cm 
thick,  just  beyond  the  outer  surface  of  the  KC-135  fuselage  is  a  major  cause  of  airborne  image 
degradation.  The  measurements  were  made  using  an  interferometer,  which  was  originally  de¬ 
signed  by  D.  Kelsall  for  another  purpose,  and  which  fortunately  became  available  in  time  for  its 
inclusion  in  the  planning  of  a  boundary-layer  experiment.  The  description  of  the  experimental 
setup  and  the  salient  results,  some  of  which  are  used  here,  are  discussed  in  Sec.  II-C  of  this 
report.  Here,  we  compare  one  of  these  results  with  the  residual  B-term  (Ref.  1.,  p.  42)  in  the 
experimental  formula  for  the  90-percent  image  diameter;  this  term  was  tentatively  ascribed  to 
boundary- layer  turbulence.  Experimentally,  the  magnitude  of  B  is  approximately  18  X  10-f>rad. 
The  point  spread  function  of  the  boundary  layer  can  be  found  by  extracting  the  Fourier  transform 
of  the  single-pass  boundary- layer  MTF  measured  at  the  LFLC  window  (see  Sec.  II-C).  Figure  11-13 
illustrates  the  MTF  as  far  as  the  diffraction  limit  of  the  1-inch  probe  beam  will  allow  its  compu¬ 
tation.  The  MTF  was  computed  by  dividing  the  diffraction-limited  MTF  into  the  "measured"  MTF 


Fig.  11-13.  Single-poss  boundary  loyer  MTF  with  diffraction  limit  divided  out. 
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(slightly  corrected  for  noise  at  the  high  spatial  frequencies),  which  included  both  boundary  layer 
and  diffraction.  Unfortunately,  the  MTF  rapidly  becomes  indeterminate  near  the  diffraction 
limit,  as  indicated  by  the  error  bars  which  are  rough  estimates  of  the  uncertainties.  The  curve 
was  extended  to  high-frequency  cutoff  by  an  arbitrary  linear  extrapolation,  which  is  probably  not 
unreasonable,  but  obviously  not  very  accurate. 

The  transform  was  carried  out  numerically  using  the  Cooly-Tukey  fast  Fourier  transform  al¬ 
gorithm.  The  resultant  point  spread  function  (Fourier  amplitude  squared)  is  shown  in  Fig.  11-14 
a.ong  with  the  90  percent  width,  which  indicates  the  points  on  a  2-dimensional  curve  which  corre¬ 
spond  to  the  circle  on  a  3-dimensional  curve  of  rotation  which  enclosed  90  percent  of  the  volume. 


Fig.  11-14.  Boundary  layer  point  spread  function  with  indication  of  90-percent  image  width. 

The  abscissa  can  be  converted  to  an  angular  scale  by  dividing  by  the  focal  length  (8*  mm)  of 
the  entrance  optics  to  the  interferometer.  The  90-percent  angular  width  was  computed  to  be 
17.5  X  10  rad,  which  is  close  enough  to  18  X  10"6  to  leave  little  doubt  that  the  B-term  repre¬ 
sents  the  boundary- layer  component. 

A  number  of  persons  participated  in  various  aspects  of  these  experiments;  among  them  are 

D.  G.  Kocher,  D.  Kelsall,  E.  A.  Davidson  and  Lt.  R.  Stiles. 

B.  W.  Bryant 

E.  TUNABLE  INFRARED  LASER  SPECTROSCOPY  OF  SULFUR  DIOXIDE 
1.  Introduction 

Sulfur  dioxide  (S02 )  is  one  of  the  most  common  and  harmful  of  the  atmospheric  pollutant 
gases.  With  the  advent  of  laser  techniques  for  remote  or  in  situ  monitoring  of  S02  and  other 
gases  by  their  1R  signatures,"  knowledge  of  their  fundamental  spectra  becomes  imperative. 

As  we  described  earlier,7  tunable  semiconductor  diode  lasers  are  being  developed  for  air 
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pollution  monitoring  applications  as  well  as  for  ultra- high- resolution  IR  spectroscopy.  Doppler- 
limited  spectra  have  already  been  obtained  for  several  gases  in  different  spectral  regions:  sulfur 
hexafluoride  (SF^),  ammonia  (NHj),  and  ethylene  (C2H4)  around  10.6  pm  using  PbQ  ggSnQ  12Te 

tunable  diode  lasers;8  and  carbon  monoxide  (CO),  nitric  oxide  (NO),  and  water  vapor  in  the 

9 

4- to  6-pm  region  using  PbS.  Se  diode  lasers. 

j  n  1  "X  X 

In  an  earlier  report  we  showed  preliminary  diode-laser  studies  of  the  band  of  SC>2, 
centered  around  8.7  pm.  The  purpose  of  the  present  report  is  to  show  agreement  between  the 
experimental  data  and  theoretical  predictions  based  on  accurate  microwave  absorption  data;  and 
thereby  to  establish  a  new,  more  accurate  value  for  the  band  center.  We  also  show  how  these 
results  can  be  used  to  predict  the  IR  spectra  of  this  gas  under  varying  conditions  of  temperature 
and  pressure,  so  that  optimum  wavelengths  for  monitoring  applications  can  be  selected. 

2.  Diode  Laser  Fabrication  and  Experimental  Procedure 

In  order  to  produce  diode  lasers  which  emit  between  1100  and  1200cm  *,  single  crystals  of 
Pb^  xSnxTe  w°re  grown  from  the  vapor  phase  by  a  closed-tube  process/1  at  three  different 
compositions:  x  =  0,0575,  0.0700  and  0.0745.  The  crystals  were  cleaved  into  rectangular  paral¬ 
lelepipeds,  each  containing  a  p-n  junction,  of  average  size  0.084  x  0.025  X  0.020  cm  for  the  --ix 
diode  lasers  used  in  this  study.  Low-resistance  ohmic  contacts  were  formed,  and  the  lasers 
were  mounted  onto  the  cold-finger  of  a  liquid-He  Dewar.  Fine-tuning  of  the  wavelength  within  a 
single  longitudinal  mode  of  the  laser  was  achieved  by  changing  the  diode  current,  which  alters 
the  refractive  index  of  the  semiconductor  through  heating  of  the  junction  region. 

Each  of  the  six  lasers  exhibited  a  maximum  total  CW  power  from  the  front  face  in  the 
6  to  44- pW  range.  Absorption  spectroscopy  was  performed  by  transmitting  the  tunable  laser 
radiation  through  a  30-cm-long  cell  containing  the  SC>2,  through  a  grating  spectrometer  and  onto 
a  Ge:  Cu  detector.  (The  purpose  of  the  spectrometer  is  to  provide  approximate  calibration  of 
the  laser  wavelength  and  to  limit  the  IR  bandpass  to  a  single  longitudinal  mode. )  The  spectra 
were  displayed  on  an  X- Y  recorder,  whose  X-axis  was  driven  by  a  signal  proportional  to  the 
diode  current  (usually  quite  linear  with  respect  to  the  wavenumber  of  the  emitted  radiation).  The 
Y-axis  was  driven  by  a  signal  proportional  to  the  transmitted  IR  power,  as  monitored  on  the  de¬ 
tector  by  a  lock-in  amplifier  synchronized  to  a  mechanical  chopper  at  the  spectrometer  entrance 
slit. 

The  S02  was  obtained  from  the  Mathe.son  Company,  Inc.  of  East  Rutherford,  New  Jersey, 
with  a  quoted  purity  of  99.98%.  Gas  pressure  was  measured  with  a  capacitance  manometer  at  and 
below  1  torr,  and  with  a  Wallace-Tiernan  gauge  at  higher  pressures.  The  usual  test  sequence 
consisted  of  scans  of  S02  at  0-,  0. 1-,  1-  and  10-torr  pressure,  followed  by  NH^  for  calibration 

purposes.  Use  was  made  of  recent  precise  measurements  of  the  NH,-line  positions  in  this  re- 

12  J 

gion  made  available  to  us  by  Rao,  Curtis  and  Yin. 

3 .  Theory 

Because  of  the  thousands  of  relatively  strong  rotational  transitions  in  SC>2,  identification  of 

individual  lines  in  the  v.  vibration- rotation  band  is  somewhat  difficult.  Fortunately,  there  are 

1  13 

accurate  microwave  absorption  data  for  both  the  ground  and  the  first  excited  vibrational  states. 

For  both  vibrational  states,  the  microwave  measurements  can  be  used  to  evaluate  the  constants 
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occurring  in  a  symmetrized  expansion  of  the  rotational  Hamiltonian  in  powers  of  the  angular 
momenta  J.14  For  values  of  J  up  to  60,  the  energy  levels  as  well  as  the  transition  frequencies 
and  strengths  for  the  Q_3  y  Q3  _y  y  P_1  _y  R_1  3  and  _3  transitions  have  been  calcu¬ 
lated.  (P,  Q  and  R  have  their  usual  meaning,  with  the  first  and  sec  tnd  subscripts  referring  to 
AK.  and  AKP,  respectively.)  In  order  to  calculate  band  strengths  from  measurements  on  a 
single  spectral  line  at  room  temperature,  the  rotational  partition  function  was  evaluated  and 
found  to  be  =  5805. 

4.  Results 

Figure  11-15  shows  a  1.  4-cm  *-wide  scan,  between  1156.3  and  1157.7cm  1,  for  SO^  at  pres¬ 
sures  of  0.1,  1  and  10  torr,  together  with  two  NH^  lines  for  calibration.  The  1  0-torr  curve 
shows  definite  pressure-broadening  and  overlap  of  closely- spaced  lines  separated  at  the  lower 
pressures.  Absorption  at  0.1  torr  is  difficult  to  observe  because  of  the  noise  level.  However, 
it  can  be  seen  for  several  of  the  strong  transitions.  Some  spectra  in  other  wavelength  regions 
where  the  laser  power  was  greater  were  considerably  less  noisy,  and  stronger  transitions  than 
those  shown  in  Fig.  11-15  could  be  easily  seen  at  pressures  below  0.1  torr.  The  individual  SO., 
transitions  were  identified  by  comparing  their  relative  positions  and  intensities  with  theory,  and 
some  of  the  stronger  lines  in  Fig.  11-15  are  tabulated  in  Table  II- 1.  Less  than  5  percent  of  the 
lines  strong  enough  to  be  measured  by  the  20  separate  spectral  regions  covered  with  the  tunable 
diode  lasers  could  not  be  explained  by  theory.  In  order  of  decreasing  probability,  we  estimate 
these  lines  to  be  caused  by:  (a)  "hot"  bands,  the  vz  state  being  8-percent  occupied  at  room  tem¬ 
perature,  (b)  other  isotopes  and  (c)  weaker  transitions  not  included  in  the  calculation. 

1 1  > !  1 1 1  , mu  ill  fill 


Fig.  11-15.  Diode  laser  scans  of  SC>2  at  pressures  of  0,  0.1,  1  and  10  torr  at  room  temperature 
with  a  30-cm  cell .  Scan  of  NHj  ot  1  torr  was  olso  made  for  calibration.  Diode  loser  excita¬ 
tion  (ond  tuning)  current  is  shown  by  upper  horizontal  scale.  Wavenumber  scale  represents 
theoretical  positions  of  observed  lines  (see  Toble  1 1— 1 )  using  1151.71  cm"'  as  band  center. 
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Fig.  11-16.  Comporison  between  measured  and  predicted  strengths  and  positions  for  some  lines  in  ^  bond 
of  SOo.  Scons  were  computer-generated  using  Lorentzian  line  profiles.  Lines  designoted  as  a,  b  ond  c 
represent  transitions  28 7  21  —  29^  24/  348,26  "*  "^7,29  ond  273,25  —  264, 22<  respectively,  which  op- 
pear  to  be  shifted  slightly  to  lower  wovenumber  (see  text).  Line  d  is  not  predicted  by  theory. 
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^  The  present  results  yield  a  value  of  1151.71  ±  0.01  cm'1  for  the  i^-band  center  frequency 
of  S  02,  which  can  be  compared  with  1151.38cm'1  obtained  by  Shelton,  Nielsen  and  Fletcher.15 
The  uncertainty  of  ±0.01  cm  1  in  our  determination  is  due  in  part  to  possible  errors  in  the  NH^ 
spectra,  the  theoretical  calculations,  and  the  intercomparison  between  the  NH^  and  SC>2  traces. 

Figure  II- 16  represents  a  detailed  comparison  between  experimental  data  from  a  diode  laser 
scan  in  the  1147-cm  region  and  the  theoretical  prediction  for  lines  in  this  region.  The  curves 
marked  "experimental"  are  Lorentzian  fits  to  the  line  intensities,  incorporating  a  linewidth  which 
best  approximates  the  shape  of  the  experimental  curve  for  the  1-torr  pressure  at  which  the  data 
were  obtained.  (Actually,  at  1-torr  pressure  the  S02  line  profile  is  a  Voigt  convolution  of 
Lorentzian  and  Gaussian  shapes.)  There  are  minor,  but  significant  differences  between  the  two 
curves,  as  indicated  at  a,  b  and  c,  where  transitions  of  high  J  value  appear  to  be  shifted  to 
the  left  in  the  experimental  curve  by  a  few  MHz/unit  change  in  J.  At  point  d  is  a  line  which  is 
not  predicted  theoretically  and  may  be  caused  by  one  of  the  reasons  mentioned  above. 

Figure  11-17  illustrates  the  pressure-broadening  rates  for  two  SC>2  lines.  Self-broadening 
of  the  153  13  -*  154  12  line  at  1163.297cm  1  occurs  at  a  rate  of  40  ±  5  MHz/torr  for  the  full 
width  at  half- maximum  at  pressures  above  a  few  torr.  Similarly,  the  effect  of  increasing  the 
partial  pressure  of  air  on  an  S02  line  at  fixed  pressure  is  shown  for  the  8i  ^  -*  8Q  transition  at 
1148.810cm  ,  where  the  broadening  rate  is  only  12  ±  2  MHz/torr.  Although  Fig.’ll-17  does  not 

represent  linewidth  measurements  on  the  same  transition,  measurements  of  the  widths  of  sev¬ 
eral  lines  do  not  indicate  any  noticeable  variation  from  the  pressure-broadening  rates  shown 
here.  Consequently,  the  ratio  of  3.  3  ±  0.7  between  self-broadening  and  air-broadening  rates  is 
significant  and  can  be  compared  to  the  5  ±  1  obtained  by  Burch,  et  al.,  when  dry  N2  is  used  in¬ 
stead  of  air.  (The  presence  of  water  vapor  in  the  laboratory  atmosphere  used  for  our  air- 
broadening  studies  may,  in  fact,  be  the  reason  for  the  smaller  ratio. ) 

By  using  measured  absorption  constants  and  pressure-broadened  linewidths  at  10-torr  pres¬ 
sure  for  several  S02  lines,  we  have  calculated  their  fundamental  line  strengths.  On  the  basis 
of  these  deter i.T  nations,  and  using  the  theoretical  ratios  of  total  partition  function  to  line  strength, 
the  total  strength  of  the  band  is  found  to  be  358  ±  20  X  lo’20cm’1  molecule'1  cm2.  This  can 
be  compared  with  the  value  of  371  ±  20  X  lo"20cm"1  molecule"1  cm2  measured  by  Burch,  et  al.,16 
and  is  in  the  range  found  by  earlier  workers1  7  using  conventional  techniques.  An  effective 
dipole  moment  of  0.086  ±  0.003  Debye  is  found  by  a  similar  calculation. 

5.  Application  of  Theoretical  Spectra  to  Pollutant  Gas  Monitoring 

Knowledge  of  the  strengths,  wavelengths  and  widths  of  individual  vibration-rotation  lines  of 
gaseous  pollutants  can  be  used  to  predict  the  IR  spectra  corresponding  to  a  variety  of  pressure 
and  temperature  conditions.  For  the  case  in  which  S02  in  smokestack  effluent  is  to  be  monitored, 
the  gas  is  at  atmospheric  pressure  and  can  have  wide  temperature  fluctuations.  Using  the  theo¬ 
retical  data  for  the  thousands  of  lines  in  the  8.  7-pm  band  of  S02,  it  was  possible  with  the  aid  of 
a  digital  computer  to  display  the  absorption  spectra  for  several  temperatures.  A  segment  of 
this  band,  from  1110  to  1142  cm  i,  is  illustrated  in  Fig.  11-18  for  selected  temperatures  from 
0  to  400  C.  Below  1126  cm  ,  absorption  is  highest  for  higher  gas  temperatures,  whereas  the 
converse  is  true  above  1126  cm  1.  Consequently,  IR  absorp  on  of  SC>2  has  the  smallest  temper¬ 
ature  dependence  around  1126  cm  ,  which  makes  this  region  particularly  appropriate  for  in-stack 
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Fig.  11-18.  Computed  spectra  far  o  portion  of  t'j  band  of  SO2  in  air  ot  atmospheric  pressure, 
ot  temperatures  from  0°  ta  400°C.  A  constant  width  of  0.3  cm-'  was  used  far  all  lines. 

monitoring  of  SOz.  Moreover,  by  alternately  tuning  the  laser  from  1125.8  to  1126.8  cm-1  and 
electronically  measuring  the  logarithm  of  the  ratio  of  th<  transmitted  signals  at  the  two  wave¬ 
lengths,  the  resulting  voltage  is  directly  proportional  to  pollutant  concentration  and  is  not  affected 
by  absorption  or  scattering  processes  which  are  invariant  over  the  small  wavelength  interval. 
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F.  EXTENSION  OF  TUNABLE  LASER  TECHNIQUES 
TO  AIR  POLLUTION  DETECTION 

We  have  continued  to  apply  techniques  using  tunable  infrared  lasers  to  the  problem  of  de¬ 
tection  and  measurement  of  air  pollution.  As  part  of  this  effort  we  have  considered  the  use  of 
acousto-optic  absorption  cells  in  conjunction  with  discretely  tunable  gas  lasers  and  tunable  diode 
lasers.  We  have  used  tunable  diode  lasers  for  in-situ  measurements  of  auto  exhaust,  as  well  as 
for  long  path  studies  of  background  concentrations.  The  work  described  concerns  the  gas 
ethylene  (C^l^),  although  we  have  also  studied  SF^,  NO  and  CO. 


\utlmr  not  .it  l.iitcoln  I  ..ilmnitory. 
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1.  Spectrophone  Inves  igation 

As  an  alternate  method  to  tunable  semiconductor  laser  transmission  measurements  for 
pollutant  deteeti'  n,  we  have  considered  a  technique  whereby  an  acoustic  signal,  oriven  by  the 
absorption  of  lit  radiation,  is  used  to  measure  trace  gas  concentration.  The  technique,  first 
investigated  in  the  19th  Century,20  has  been  used  recently  to  measure  decay  times  in  gases 
and  to  detect  traces  of  methane  and  nitric  oxide  in  air  samples. 

The  acousto-optic  absorption  cell,  or  spectrophone,  can  assume  a  variety  of  lorms.  The 
constructions  we  have  used  consist  of  an  inner  and  an  outer  chamber  separated  by  a  capacitance 
microphone.  The  air  sample  is  admitted  to  the  outer  chamber  and  diffuses  into  the  inner  cham¬ 
ber  by  means  of  a  small  hole.  Laser  radiation  passes  in  and  out  of  windows  attached  to  the 
inner  chamber.  If  the  laser  frequency  coincides  with  an  absorption  line  of  the  gas,  some  of  its 
energy  is  lost  to  the  gas  and  appears  as  translational  energy  causing  a  pressure  differential 
between  the  inner  and  outer  chambers.  V.y  chopping  the  laser  beam  and  synchronously  detecting 
the  capacitance  change  produced  by  this  pressure  differential,  we  have  measured  ethylene  con¬ 
centrations  as  low  as  28ppb.  This  measurement  used  a  0.5  W  stable  C02  laser  tuned  to  the  P(14) 
line,  which  nearly  coincides  with  an  absorption  line  in  ethylene. 

Unfortunately,  at  this  low  concentration  a  signal  produced  by  absorption  in  the  windows  is 
dominant  and  must  be  subtracted  from  the  measured  signal  to  calculate  the  effect  of  the  trace 
gas.  In  this  particular  example,  the  (BaF2)  window  signal  was  2.5X  the  intensity  of  the  signal 
due  to  28  ppb  of  C2II4.  Consequently,  a  high  degree  of  stability  in  laser  amplitude  and  in  con¬ 
figuration  alignment  is  required  to  complete  the  measurement. 

This  is  not  the  case  with  lower  power  lasers  where  the  sensitivity  of  the  spectrophone  deter¬ 
mines  the  lower  limit  of  detection.  Several  of  the  spectrophone  cells  we  constructed  were 
capable  of  detecting  less  than  lo'8W  of  absorbed  power  (in  the  absence  of  a  large  window  absorp¬ 
tion  signal).  With  a  5-pW  PbSnTe  diode  laser  we  were  thus  able  to  detect  an  a'el.  of  2  x  10 
where  the  length  L  used  was  10cm,  c  is  the  gas  concentration  (ppm)  and  o'  is  its  absorption 
coefficient  in  (cm  -  ppm)'1.  Since  a'  is  typically  1  o'6 /cm  -  ppm  for  0,1^,  this  corresponds  to 
a  minimum  detectable  concentration  near  200  ppm. 

In  Fig.  11-19  we  show  a  comparison  between  the  spectrophone  output  and  the  signal  received 
by  a  Ge:  Cu  detector  after  the  5-p.W  PbSnTe  laser  radiation  has- passed  through  the  spectrophone. 
The  diode  radiation  was  tuned  over  a  region  near  10,6  pm,  and  the  Doppler-limited  spectrum  is 
evident  in  both  the  spectrophone  and  the  transmission  signals.  While  the  minimum  detectable 
concentration  using  the  spectrophone  was  more  than  102ppm,  a  detection  limit  ot  1  ppm  was 
possible  with  the  liquid-Ue-cooled  Ge:Cu  detector  under  these  conditions.  This  reflects  the  fact 
that  the  spectrophone  NEP  was  near  1  0  8  W,  while  the  GeCu  detector  had  a  NEP  of  approximately 
3  X  1  O'11  W. 

2.  In -Situ  Auto  Exhaust  Measurements 

The  diode  laser,  being  compact,  stable  and  easily-operated,  is  highly  suitable  for  field 
measurements.  Utilizing  this  capability  we  have  performed  fast  response  in-situ  measurements 
of  the  ethylene  concentration  in  auto  exhaust.  The  undiluted  exhaust  was  drawn  through  a  win¬ 
dowless  1.1 -m  sample  cell  and  exhausted  with  a  standard  garage  exhaust  fan.  Radiation  from  a 
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Fig.  11-19.  Detection  of  C2H4  as  diode  laser  emission  wovelength  is  tuned  neor  10.6  pm: 
(o)  transmission  through  spectrophone  cell  (detection  t  =  0. 03  sec);  (b)  spectrophone output 
(detection  t  =  3  sec).  Cell  length  =  10  cm;  pressure  =  5  torr;  oir  pressure  =  15  torr. 


Fig.  11-20.  On-line  CoH^  emission  test  for  1972  station  wogon  (system  using  pulsed 
Pbg  ggSng  diade  loser  at  10.57  pm). 
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PbSnTe  diode  laser  operating  at  10.57  pm  was  collimated,  passed  through  the  sample  cell,  and 
collected  by  a  Ge:Cu  detector.  The  detector  output  was  amplified  with  a  logarithmic  amplifier 
and  fed  to  a  strip  recorder  whose  deflection  was  proportional  to  concentration. 

Calibration  was  accomplished  by  substituting  a  cell  with  known  ethylene  concentration.  Cali¬ 
bration  and  lack  of  interference  from  other  gases  was  further  checked  by  low  pressure  (Doppler 
regime)  studies  of  exhaust  samples  and  samples  of  known  concentration. 

A  typical  measurement  sequence  is  shown  in  Fig.  11-20.  In  Fig.  II-20(a)  the  C2H4  emission 
from  a  19  72  station  wagon  is  measured  after  a  "normal"  start.  The  emission  shows  two  con¬ 
centration  peaks  after  which  the  level  approaches  350  ppm.  Figure  II- 2 0(b )  shows  a  "rich" 
start  for  which  the  accelerator  pedal  was  depressed  once  before  starting.  The  emission  exceeds 
2000ppm  briefly  then  approaches  the  350  ppm  steady- state  level.  Removal  of  the  air  cleaner 
on  this  vehicle  leaned  the  fuel-air  mixture  and  reduced  the  steady  state  ethylene  concentration 
to  l60ppm.  Other  vehicles  tested  showed  somewhat  higher  concentrations.  Figure  11-21  shows 
the  ethylene  content  of  a  poorly  tuned  1964  sedan.  The  average  concentration  is  800 ppm.  The 
wide  trace  is  not  due  to  system  noise,  but  rather  to  rapid  fluctuations  in  this  car's  ethylene 
output. 

These  tests  were  designed  to  demonstrate  the  suitability  of  diode  laser  transmission  studies 
to  real-time  measurement  of  pollution  sources.  They  will  be  extended  to  include  other  gases, 
such  as  nitric  oxide  and  carbon  monoxide. 


fifY* 


Fig.  11-21.  On-line  C2H4  emission  test  for  untuned  1964  sedan  (system  using  pulsed 
Pbg  ggSn^  ^Te  diode  laser  ot  10:57  pm). 
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Fig.  11-23.  Measurement  of  ambient  C2H4  concentration  using  Pbg  ggSng  ]2  Te  diode 
laser  (0. 5-km  path). 
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3.  Long-Path  Laser  Transmission  Studies 

To  increase  the  sensitivity  of  the  d -ode  laser  transmission  technique  to  measure  pollutant 
levels  in  ambient  air  (typically  10  to  100  ppb)  it  is  necessary  to  use  longer  path  lengths  than  are 
convenient  for  laboratory  application.  Also,  the  integrated  concentration  obtained  from  long 
path  measurements  is  often  more  suitable  than  point  by  point  results  for  regional  air  quality 

studies. 

We  have  operated  the  PbSnTe  diode  laser,  described  above,  over  a  0.5-km  path  from  Lincoln 
Laboratory  to  a  remote  reflector  and  back.  The  experimental  path,  illustrated  in  Fig.  11-22, 
passed  over  the  Laboratory  parking  lot  and  allowed  us  to  measure  changes  in  the  ethylene  content 
during  peak  traffic  hours.  Figure  11-23  shows  an  increase  in  the  measured  C2H4  concentration 
of  over  1  ppm  during  the  late  afternoon  departure  of  Laboratory  personnel.  Peaks  due  to  individ¬ 
ual  cars  (such  as  the  poorly  tuned  1964  sedan  mentioned  above)  were  also  noted. 

In  this  experiment,  the  laser  radiation  was  collimated  to  a  2  x  4-cm  beam  before  transmis¬ 
sion.  It  was  collected  by  a  larger  mirror,  used  off-axis,  and  concentrated  on  the  liquid- He- 
cooled  detector.  With  this  very  simple  arrangement,  we  were  able  to  collect  more  than  30  per¬ 
cent  of  the  transmitted  energy.  The  SNR  obtained  was  limited  by  atmospheric  fluctuations, 
rather  than  diode  power  or  detector  considerations.  As  a  result  of  these  measurements,  we 
expect  to  be  able  to  extend  the  path  length  to  many  kilometers. without  difficulty.  By  using  the 
two-pulse  frequency  scanning  technique,  developed  for  the  EPA  smokestack  measurement  appara¬ 
tus,  we  expect  to  obtain  detection  sensitivities  near  1  ppb. 

H.  A.  Pike 
E.  D.  Hinkley 
A.  R.  Calawa 
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III.  LASER  RADAR  AND  TRACKING 


A  INSTRUMENT  FOR  ATMOSPHERIC  MODULATION  TRANSIT  £R 

FUNCTION  MEASUREMENTS 

1,  Introduction 

A  modulation  transfer  funetion  (MTF)  instrument  which  is  of  simple  design  and  easily 
transported  has  been  constructed  for  visible  atmospheric  propagation  measurements.  A  picture 
of  this  instrument  and  a  schematic  drawing  of  its  operation  are  shown  in  Figs.  III-Ka.b). 

An  8.6-cm-diameter  Questar  telescope  is  used  to  form  an  image  of  a  distant  point  source, 
e  g  a  He-Ne  laser.  This  point  source  is  re-imaged  onto  the  rotating  MTF  reticle,  which  is  a 
Siemann  star  pattern  of  alternately  reflecting  and  nonreflecting  2“  segments  or  bars.  The  bar 
spacing,  that  is,  the  spatial  frequency,  is  a  function  of  radius  on  the  reticle  and  is  selected  by 
tilting  the  scan  mirror.  The  modulation  of  the  refleeted  light  at  each  radius  setting  produced  by 
the  spinning  reticle  and  detected  at  the  photomultiplier  gives  a  direct  measure  of  the  MTF  at 

that  spatial  frequency. 


i|-5-40»1 


-5- 4039  I 
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p;g. IH-1.  MTF  instrument:  (a)  photograph;  (b)  schematic  diagram. 
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Fig.  1 11-2.  Instrumental  modjloiion 
transfer  function. 


Figure  III-2  shows  the  theoretically  calculated  MTF  for  the  instrument  along  with  a  corre¬ 
sponding  experimental  measurement.  The  particular  shape  is  due  to  the  central  obstruction  in 
the  Questar  aperture  (tj  =  vjvz  =  0.39).  Also  shown  in  Fig.  III-2,  is  an  example  of  the  instru¬ 
mental  MTF  degraded  by  the  atmosphere.  The  case  shi  wn  is  for  a  5.4-km  horizontal  path  and 
was  taken  at  13:10  EDT.  The  results  of  a  more  comprehensive  set  of  measurements  on  the 
atmospheric  MTF,  made  in  conjunction  with  10.6-gm  propagation  measurements  also  described 
herein,  will  be  presented  in  the  next  Optics  Research  Report. 

J.  Z.  Holtz 
T.  J.  Gilmartin 
L.  J.  Prugnarola 

2.  MTF  Electronics 

The  electronics  used  with  the  MTF  instrument  consists  of  a  square-wave  generator  and  a 
signal  processor. 

The  square-wave  generator  consists  of  a  light  source  and  photodetector  in  an  integral  hous¬ 
ing  with  a  beam  splitter  and  lens  such  that  the  light  is  focused  on  the  rotating  reticle  and  the  re¬ 
turn  beam  is  focused  on  the  photodetector.  Following  the  photofet  detector  are  an  amplifier  and 
an  emitter  follower  which  provides  a  low  output  impedance.  The  square-wave  generator  pro¬ 
vides  the  reference  signal  when  a  lock-in  amplifier  is  used  to  process  the  data. 

The  signal  processor  provides  an  output  proportional  to  the  percentage  modulation  im¬ 
pressed  on  the  incoming  beam  by  the  rotating  reticle.  It  includes: 

(a)  Photomultiplier 

(b)  Photomultiplier  buffer 

(c)  Low-pass  active  filter  before  the  signal  strength  meter 

(d)  AGC  loop  amplifier  and  filter,  which  control  the  voltage  applied  to 
the  photomultiplier  and  thus  its  gain 
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(e)  Active  bandpass  filter  tuned  to  the  frequency  produced  by  the  bar 
pattern  of  the  rotating  reticle 

(f)  Full-wave  rectifier  and  filter,  which  convert  the  output  of  the  bandpass 
filter  into  a  DC  voltage  proportional  to  the  percentage  modulation. 


Fig.  111-3.  MTF  control  ond  signol  processing  electronics. 

Outputs  are  available  at  various  points  for  recording  or  for  connecting  other  measuring  instru¬ 
ments.  Figure  III-3  is  a  block  diagram  of  the  electronics.  M>  D  Timmerman 

B.  AUTOMATIC  RANGE  TRACKER 

Detailed  plans  for  an  automatic  range  tracker  for  the  Firepond  laser  radar  installation  were 
devised.  Implementation  of  the  plans  has  been  cancelled  because  of  planned  program  and  system 
changes  at  Firepond. 

An  automatic  switch  designed  to  interchange  the  tuning  relationships  of  the  local  oscillator 
(L.O)  to  the  master  oscillator  (MO)  as  controlled  by  the  delta  f  loop  servo  has  been  constructed 
and  is  ready  to  be  installed  and  tested.  It  will  switch  the  output  frequency  of  the  LO  from  5  Mllz 
below  the  MO  frequency  to  5  MHz  above  the  MO  frequency  when  the  Doppler  frequency  goes  to 
zero.  In  this  manner,  the  signal  frequency  at  the  detector  mixer  output  will  be  prevented  from 
going  to  zero  regardless  of  whether  the  moving  target  is  on  the  incoming  or  outgoing  leg  of  its 
pass.  Figure  III-4  is  a  block  diagram  of  the  switching  scheme. 

Design  and  construction  of  the  remote  control  circuits  to  be  used  with  the  AIL  monopulse 
tracking  receiver  have  been  completed.  The  circuits  are  ready  to  be  installed  and  tested  at 
Firepond  when  the  receiver  arrives. 


M.  D.  Zimmerman 
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Fig.  111-4.  Automotic  Doppler  switchover. 

C.  ATMOSPHERIC  OPTICAL  COHERENCE  MEASUREMENTS  AT  10.6  AND  0.63  (ira 

The  performance  of  large  aperture  optics  (1)  >  I  m)  at  10. 6  pm  in  projection  and  heterodyne 
reception  has  previously  been  estimated  by  extrapolating  from  measurements  made  at  shorter 
wavelengths.  Since  large  and  complex  systems  are  being  contemplated  at  10.6  pm,  it  is  of  great 
practical  interest  to  measure  directly  the  effect  which  atmospheric  refractive  inhomogeneities 
will  have  on  such  systems.  To  this  end,  a  48-inch-diameter  (1.2  m)  Cassegrainian  telescope 
was  focused  at  a  range  of  5.4km,  where  the  10.6-pm  spot  fell  on  a  diffuse  metal  surface,  was 
scanned  with  an  IK  imaging  camera  and  recorded  for  later  computer  analysis.  Simultaneously, 
the  optical  quality  of  parallel  paths  1  and  2  m  distant  from  the  primary  path  were  monitored  at 
0.63  pm  with  a  MTF  instrument  and  with  a  two-point,  white  light  source  resolution  instrument. 

The  optical  paths  were  horizontal  with  an  average  altitude  of  68  m. 

Data  have  now  been  taken  under  a  variety  of  meteorological  conditions  and  over  several 
diurnal  cycles.  However,  while  most  of  the  data  from  the  visible  instruments  have  been  analyzed, 
the  computer  analysis  of  the  10.6-pm  data  will  not  be  completed  for  several  weeks.  For  this 
reason,  a  detailed  description  of  the  experiment  or  of  the  results  will  not  be  presented  at  this 
time,  but  some  qualitative  observations  will  be  made. 

Defining  the  atmospheric  coherence  diameter  rQ  to  be  the  diameter  which,  when  used  in  a 
vacuum,  allows  the  same  resolution  as  docs  the  atmosphere  when  not  aperture  limited,  the  at¬ 
mospheric  r  was  measured  by  the  MTF  instrument  to  range  from  about  -j  to  20cm  at  0.63  pm. 

O  L  /c 

Since  r  is  predicted  to  scale  as  \  ,  r  (10.6  pm)  is  expected  to  equal  29.4  r  (0.63  pm)  and  the 

o  O  o 

rn(10.6pm),  implied  by  the  visible  measurements,  ranges  from  30cm  to  6  m.  Indeed,  the  image 
projected  by  the  1.2-m  aperture  at  1  0#  pm  was  observed  to  be  almost  perfectly  diffraction  lim¬ 
ited  occasionally,  near  diffraction  limited  most  often,  and  badly  degraded  occasionally.  A 
measurement-by-mcasuroment  correlation  of  coherence  data  taken  at  the  two  wavelengths  is 
presently  being  made. 

The  10.6-pm  spot  wander  was  observed  to  range  from  a  fraction  of  the  diffraction- limited 
spot  radius  to  several  spot  radii  at  extremes  under  typical  wind  and  turbulence  conditions.  The 
correlation  time  for  spot  wander  and  image  distortion  was  typically  a  fraction  of  a  second. 


88 


Section  111 


A  more  detailed  and  comprehensive  description  of  the  results  of  these  propagation  measure 


ments  will  be  given  in  the  next  Optics  Research  Report. 


T.  J.  Gilmartin 
J.  Z.  Holtz 
R.  N.  Capes 


D.  InSb  ISOLATOR 

W,  have  previously  reported  that  when  10.6-pm  radiation  is  incident  on  an  InSb  sample  that 
is  rigidly  attached  to  a  heat  sink,  beam  distortion  occurs.  The  stress  caused  by  local  heating 
will  not  be  uniform  in  this  case  and  the  resultant  unsymmetrical  change  in  the  index  of  refraction 

will  cause  the  distortion. 

A  new  holder  has  been  designed  which  allows  the  sample  to  expand  uniformly.  The  holder 
is  shown  in  Fig.  I1I-5.  Tests  u  ling  a  12-W  beam  incident  on  a  4-mm  piece  of  InSb  held  by  'he 
new  holder  show  none  of  the  distortion  previously  observed. 


Fig.  1 1 1-5.  InSb  holder. 


A  profilometcr  was  constructed  which  scans  the  beam  at  l*-mch/min. 
through  a  20  to  30  thousandths  aperture.  The  time  constant  of  the  detector 
from  the  profilometer  show  no  gross  differences  between  the  beam  directly 


(25  thousandths/sec) 
is  10'3  sec.  Results 
from  the  laser  and 


one  that  has  passed  through  an  InSb  sample. 

An  isolator  was  constructed  using  Faraday  free  carrier  rotation.  Pictures  of  the  system 
are  shown  in  Figs.  III-6  and  -7.  The  Dewar  and  a  polarizer  arc  shown  in  Figs.  III-8  and  -9, 


respectively. 

The  isolator  has  been  tested  in  the  laboratory. 
Material 
Concentration 
Thickness 

Verdet  constant 
Field  for  45  °  rotation 
Operating  temperature 


Characteristics  of  the  10.6-pm  isolator  are: 

n  —  InSb 

1.0  x  lO1^  cm"3 
4  mm 

1.1  X  10-2  deg/cm-G 

10.2kG 

77  °K 
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Fig .4*1 1—8 .  Exploded  view  of  Dewor:  on  left, 
the  outer  Dewor  with  Ge  A-R  cooted  windows; 
on  right,  the  inner  liquid  nitrogen  reservoir  with 
InSb  ottoched.  Linuid  nitrogen  flows  to  the  cop¬ 
per  plug  visible  just  below  the  InSb  holder. 


PH0-5S! 


Fig.  111-9.  Exploded  view  of  Ge  brewster  ond  polarizer. 
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Insertion  loss 
Isolation 

Tc..al  reverse  loss 

Maximum  forward 
throughput  power 

Safety  factor 


2  ±  0.  5  dB 
20  ±  0.  5dB 
22  ±  ldB 
30  W  CW 

2.9  (calculated) 


The  safety  factor  is  the  ratio  of  the  radius  of  any  aperture  to  the  beam  waist  at  the  worst  point 
in  the  optical  path.  In  this  case  it  occurs  at  the  point  where  the  pole  faces  and  the  magnet  coils 
meet. 

Isolation  seems  to  be  limited  by  non-uniformity  in  the  magnetic  field  across  the  i-inch 
aperture  in  the  pole  piece.  Maximum  isolation  based  on  the  limitation  of  the  Ge  brewster  angle 
polarizers  is  26dB.  (There  is  a  23  percent  transmission  of  perpendicular  polarization  through 
each  surface.  There  are  4  surfaces,  thus  the  isolation  is  0.234.]  Notice  we  are  6dB  below  this 
level  indicating  that  this  is  not  the  limiting  factor. 

Some  improvement  might  be  possible  by  using  a  magnet  with  larger  pole  faces.  Larger 
pole  faces  generally  reduce  variation  of  field  across  a  ^-inch-diameter  hole  in  the  pole  faces. 

An  investigation  of  the  power  handling  capability  of  the  isolator  when  300 -nsec  pulses  are 
incident  will  be  undertaken.  The  effect  of  InSb  on  the  waveform  of  the  pulse  will  be  examined. 


L.  R.  Tomasetta 
W.  E.  Bicknell 
D.  H.  Bates 


